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Abstract 



Dynamic and static properties of the classical Fukuyama-Lee-Rice model and the 
renormalization and phase diagrams of the related quantum model with phase-slips are studied. 
In the first part, the phase correlation function is calculated in the weak pinning limit by an 
one-loop renormalization group calculation and exactly in the strong pinning case. Further, the 
creep dynamics of these quasi-one-dimensional systems is studied by analytical and numerical 
approaches. 

In the second part, the phase diagrams of the quantum model is studied by an anisotropic, 
finite-temperature renormalization group calculation. 
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Summary 

Charge density waves (CDWs) are possible low energy states of quasi-one-dimensional metals, 
which appear below a certain transition temperature Tc, typically below 200K. The underlying 
mechanism for this phase-transition - called Peierls transition - is the coupling of the band 
electrons to the underlying lattice (electron-phonon coupling). At low temperatures the system 
can reduce its energy by formation of a periodic lattice distortion, which is of the order one 
percent of the interatomic distance. This deformation leads to a gap in the electronic dispersion 
and with that to a gain of electronic energy larger than the loss of elastic energy. 
CDWs can be described by a Ginzburg-Landau model as elastic, disordered systems. 

In the first part of this thesis (chapter 2) I study dynamic and static properties of the classical 
Fukuyama-Lee-Rice model: In the weak pinning limit the temperature dependency of the 
pair correlation function is calculated using a one-loop approximation of the related Burgers 
equation. This result is verified by a Monte-Carlo simulation. For strong disorder, the correlation 
function is calculated exactly, considering the order statistics of the positions of the impurities, 
at zero temperature. 

As the central problem, the collective dynamics, especially the creep-dynamics, is studied 
by analytical and numerical approaches. Contrary to high dimensional systems an analytical 
expression for the creep velocity or current is found, which is in very good agreement with 
numerical simulation results. This result extends the well known results for D > 2. 

In the second part (starting with chapter 3) the classical model is extended by considering 
quantum fluctuations and space-time vortices [phase-slips). 

With an anisotropic, finite-temperature renormalization group approach the phase diagrams 
for this model are studied. In former works only the zero temperature behavior for either ID 
disordered, electronic systems or systems with dislocation is studied. 

The full consideration of thermal fluctuations and phase-slips is therefore an essential extension 
of these works and leads to a new scenario for the well know delocalization transition at T = 0. 

partly published in 

Andreas Glatz and Mai Suan Li, Collective dynamics of one-dimensional charge density waves, 
Phys. Rev. B 64, 184301 (2001). (excerpts of chapter 2) 

Andreas Glatz and Thomas Nattermann, One-Dimensional Disordered Density Waves and 
Superfluids: The Role of Quantum Phase Slips and Thermal Fluctuations, Phys. Rev. Lett. 88, 
256401 (2002). (based on chapters 3 and 4) 



2 



Einfuhrung 



1 Einfuhrung 
1.1 Einleitung 

In dieser Arbeit werden einige neue Aspekte eindimensionaler Ladungsdichtewellen behandelt. 
Es werden statische und dynamische Eigenschaften mit Hilfe des klassischen Fukuyama-Lee- 
Rice-Modells, welches nur thermische und Unordnungsfluktuationen beriicksichtigt, und die 
Phasendiagramme eines erweiterten, phanomenologischen Modells, welches zusatzlich noch 
Quantenfluktuationen und Dislokationen beinhaltet, bei tiefen Temperaturen untersucht. 
Doch zunachst eine kurze Einfuhrung zu Ladungsdichtewellen (CDW) und warum diese 
interessant sind: 

Schon 1955 fand R.E. Peierls anhand von theoretischen Uberlegungen, dass eindimensionale 
Metalle welche an das unterliegende Gitter gekoppelt sind, aufgrund der Elektron-Phonon- 
Wechselwirkung nicht stabil sind und eine Ladungsdichtewelle ausbilden konnen. CDWs konnen 
auch in hoheren Dimensionen auftreten, werden aber vorwiegend als quasi-eindimensionales 
Phanomen beschrieben. 



CDW Materialien 

• NbSeg 

• TaSs 

• NbSa 

NbSe, 



• Nb 



Chain direction 
OSe 



K0.3M0O3 

RbosMoOs 

(TaSe4)2l 



• (TMTSF)2PFs 

• (MX4)2Y 

• TTF-TCNQ 




K03M0O3 




Abbildung 1: Kristallstruktur einiger CDW-Materialien. Diese bestehen aus schwach 
gekoppelten Molekiilketten. Die Ubergangstemperaturen Tp zur CDW liegen fur NbSea bei 
lAbK und (!) 59il' (diese Substanz bildet zwei unabhangige CDW Phasen aus) und fur 
K0.3M0O3 bei imK. 



Einige CDW Materialien sind in Abb. [J mit Kristallstruktur zusammengestellt. Aufgrund ihrer 
Struktur konnen sie als eindimensionale Systeme behandelt werden. In Abb. [21 und Ql sind 
Fotos von NbSe3 und K0.3M0O3 zu sehen. Am Beispiel von NbSe3 kann man erkennen, 
dass die Materialien als sehr feine Faden (whiskers) vorliegen konnen und somit deren 
transversale Ausdehnung, wenn diese eine gewisse Ausdehnung (k, Lc, s. Kap.|2I) unterschreitet, 



1.1 Einleitung 
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Abbildung2: NbSea Kristallfaden Abbildung 3: K0.3M0O3 Kristalle. (Aufgrund der 

tiefblauen Farbe auch Blau-Bronze genannt.) 

vernachlassigt warden kann. 

Betrachten wir zunachst den ID Elektron-Phonon Hamiltonoperator (Frohlich Hamiltonian) 
fur spinlose Fermionen (In Gru94| R^riiSSI ist eine ausfuhrliche Darstellung der theoretischen 
und experimentellen Grundlagen zu Ladungsdichtewellen zu finden.) 

k q k,q 

wobei c^^'' und Uq '^ Vernichtungs- (Erzeugungs-) Operatoren fur Elektronen bzw. Phononen sind 
und Ck = dq ist die Kopplungskonstante der Elektron-Phonon-WW. Mit der Definition 
des komplexen Ordnungsparameters 

A = \A\e^^^g2k,{b2k,+bl^) (1.1) 

wobei kp der Fermi-Wellenvektor ist, und der Gitter-Verschiebung u{x) = YlqQq^^'^^ '^^^ 

/ ^ \l/2 . 

Qi ~ \ 2Moj ) (^g + ^-g) ist die lonenmasse) findet man fur {u{x)) folgenden Ausdruck 

Das bedeutet, dass das Gitter, welches ohne Elektron-Phonon-WW die Gitterkonstante a hat, 
aufgrund dieser WW periodisch gestort wird mit Periode Aq = 2iT/{2kp), was zur Ausbildung 
einer Liicke (A) im Elektronenspektrum bei zbfci? fuhrt (s. Abb.|^. 
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Einfuhrung 




Abbildung4: Dispersion in einem ungestorten eindimensionalen Metall mit halber BandfLillung 
(a) und mit Peierls-Storung (b). 



Die elektronische Energie nimmt fur kleine Verschiebungen oc u^lnn ab |R573| . die elastische 
Energie mit oc zu. Insgesamt erhalt man durch die Ausbildung dieser periodischen 
Gitterstorung einen Energiegewinn. Aus der Minimalbedingung fur die Energie kann man den 
Wert von |A| bestimmen. Mit dieser Gittermodulation wird auch die Ladungsdichte mit der 
Periode Aq moduliert, wobei pro Periode zwei Elektronen vorhanden sind. 
Fiir das phanomenologische Modell, welches im nachsten Kapitel behandelt wird, entwickelt 
man einen Landau-Ginzburg-Hamiltonian fiir den Ordnungsparameter A(a;), welcher zunachst 
Fluktuationen der Amplitude |A| und der Phase (/? enthalt. Da die Fluktuationen der Amplitude 
bei tiefen Temperaturen keine Rolle spielen^, werden diese im Rahmen dieser Arbeit nicht weiter 
(direkt) behandelt. In Kapitel [31 wird unser Modell erweitert und Dislokationen in der Phase 
zugelassen. Da das Phasenfeld an den Orten der Dislokationen einen Sprung hat, muss der 
Ordnungsparameter an diesen Stellen verschwinden und damit auch die Amplitude. 
Als einzig relevanter Beitrag bleibt noch der elastische Energieterm 

We/ = y dxU^d^^f (1.2) 

ubrig (die Zeitabhangigkeit des Ordnungsparameters bzw. ein daraus folgender kinetischer 

Energieterm wurde vernachlassigt, siehe dazu Kapitel I^UHJ. Folgt man den Berechnungen in 

| Gru 94| findet man c = hvp/i'^T^)- 

^Die Amplitudenfluktuationen konnen vernachlassigt werden, da diese im Landau-Ginzburg-Hamiltonian als 
quadratische Terme auftreten, und somit bei tiefen Temperaturen unterdriickt werden. 



1.2 Bekannte Ergebnisse 
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1.2 Bekannte Ergebnisse 

In diesem Abschnitt mochte ich kurz darstellen, was im Hinblick auf die in dieser Arbeit 
behandelten Themen an Vorarbeiten geleistet wurde. 

Viele Arbeiten, die die statischen Eigenschaften des klassischen Modells behandein, sind 
schon Ende der 1970er bis Anfang der 1980er Jahre geschrieben worden. Hervorzuheben 
ist unter anderen die Arbeit von M.V. Feigel'man |Fei80 |. in der Eigenschaften von 
klassischen, ungeordenten Systemen bei tiefen Temperaturen, wie z.B. Ladungsdichtewellen 
mit Unordnungspotential oder auch Spinketten mit zufalligen Anisotropien, untersucht werden. 
In dieser Arbeit wird, neben der Berechnung der kollektiven Leitfahigkeit von CDWs, auch die 
Phasen-Paarkorrelationsfunktion im Fall von schwacher Unordnung bei T = mittels einer 
Transfermatrixberechnung bestimmt. 

J. Villain und J.F. Fernandez |VF84| haben diese Korrelationsfunktion fur T — > fur ein 
harmonisches System mit Unordnungsfeld in hoheren Dimensionen (2 < d < 4) bestimmt, 
aber auch in d = 1 im Fall eines starken Unordnungspotentials. 

Wir werden diese Ergebnisse auf endliche Temperaturen fur CDWs erweitern, und im Fall 
von starkem Unordnungspotential bei T = ein exaktes Ergebnis unter Beriicksichtigung der 
geordneten Statistik der Positionen der Verunreinigungen prasentieren. 

Die Dynamik von elastischen Medien mit Unordnung bei einer extern angelegten Kraft bzw. 
eines elektrischen Feldes ist ebenfalls schon in mehreren Arbeiten untersucht worden. 

Wichtige Arbeiten in diesem Zusammenhang in Bezug auf Ladungsdichtewellen sind u.a. 
von H. Fukuyama u nd P.A. Lee \^y^ und P.A. Lee und T.M. Rice | LK79| geschrieben 
worden, in denen der Depinning-Ubergang bei T = untersucht und eine Abschatzung 
des Depinningfeldes angegeben wird, aber auch Leitfahigkeiten mit Hilfe der Kubo Formel 
berechnet wurden. 

Eine Sammlung von neueren Arbeiten zur Dynamik, auch experimentelle, ist in [BM993 zu 
finden. 

Der Beitrag dieser Arbeit zur Dynamik beschaftigt sich im wesentlichen mit der Kriechdynamik 
eindimensionaler CDWs bei tiefen, aber endlichen Temperaturen und kleinen angelegten 
Feldern. In hoheren Dimensionen ist diese Dynamik wohl bekannt [NSOO] (fur d = 2 auch 
|T592p und wir werden dies in Abschnitt 12.31 fiir d = 1 studieren. 

In einer Dimension erhalten wir fur die Kriechgeschwindigkeit einen analytischen Ausdruck, 
wobei die Feldabhangigkeit durch eine hyperbolische Sinus-Funktion gegeben ist. Dieses 
Verhalten haben wir ausfLihrlich numerisch untersucht. Bei T = haben wir auBerdem den 
Depinning-Ubergang untersucht und den kritischen Exponenten der CDW-Geschwindigkeit 
bestimmt. Unser Ergebnis stimmt mit den von anderen Autoren |5L90I IM593| . ebemfalls 
numerisch, bestimmten Werten (innerhalb der Fehlergrenzen) gut uberein. 

Relevante Arbeiten die bzgl. der untersuchten Phasendiagramme bei tiefen Temperaturen des 
erweiterten Modells schon geschrieben wurden, werden in Kapitel |H1 mit unseren Ergebnissen 
verglichen. 
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Einfuhrung 



1.3 Uberblick 

In dieser Arbeit wird der Einfluss von Unordnung in Ladungsdichtewellen eine zentrale Rolle 
spielen, welcher in einer Dimension besonders relevant ist und die Eigenschaften des Systems 
stark bestimmt. Die Arbeit ist wie folgt aufgebaut: 

• In Kapitel [2 wird zunachst der Einfluss der Unordnung im klassischen Fall (Fukuyama- 
Lee-Rice-Modell) auf die Phasen-Paarkorrelationsfunktion {{<f{x) — Lp{x')Y) im Fall 
von schwachem und von starkem Unordnungspotential behandelt. AnschlieBend wird 
die Dynamik des Systems bei einer von auBen angelegten treibenden Kraft bzw. einem 
elektrischen Feld untersucht. 

• In Kapitel Qlwird dieses Modell um einen kinetischen Term, da bei tiefen Temperaturen 
auch Quantenfluktuationen beriicksichtigt werden mussen, und einen Term, welcher den 
Effekt von Dislokationen bzw. "phase-slips" berucksichtigt, erweitert. 

• In Kapitel |3 werden die Phasendiagramme dieses erweiterten Modells ausfuhrlich 
besprochen und 

• in Kapitel El deren physikalische Konsequenzen diskutiert und mit verwandten Systemen 
verglichen. 
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2 Klassisches Modell 

In diesem Kapitel behandele ich das klassische Fukuyama-Lee-Rice Modell fFuk76l ILR79 | 
fiir eindimensionale Ladunsdichtewellen. Zunachst wird der EinfluB der Unordung auf die 
statischen Eigenschaften, insbesondere die Paarkorrelationsfunktion {{(p{x) — ip(O)y) bei 
niedrigen Temperaturen untersucht. Dabei wird zwischen schwachem (Abschnitt I2.2.1j) und 
starkem (Abschnitt I2.2.2j) Unordnungspotential unterschieden. Mit Hilfe einer Monte-Carlo 
Simulation werden die Resultate in Abschnitt 12. 2. 31 chlieBlich uberpruft. 
In Abschnitt l2.3l wird das dynamische Verhalten des Systems bei einer von auBen an das System 
angelegten Kraft bzw. elektrischen Feld analytisch und numerisch untersucht. 



2.1 Ladungsdichte 

Als fundamentale Observable wird ein Ausdruck fiir die Ladungsdichte benotigt. 
Dazu wird einfach die Teilchendichte als p{x) = n ^ • 5(x — Xj) geschrieben, wobei Xj = 
Xi — die Positionen der Elektronen sind und kleine Abweichungen von der mittleren 
Position Xi = ia beschreibt [ifi hat demnach die Funktion eines Verschiebungsfeldes). Dies 
laBt sich wie folgt umformen: 



p{x) = n 



J dX6{X - Xi)6 (^x-X + 



n I dX \ ' -e'^^^b ( X - X + 



PSF " / ~" ^ a" " V " ' 2tt 

m 



mit if{xi) = ifi und Qm = li^mja. Fuhrt man noch die neue Variable y = X — ^f^a ein, 



I uiiiL iiiaii iiw^ii uic ii^u^ vaiiaui^ y — j\. — 

folgt (vgl. mit Herleitung in [HalSl p 

a \ zir / — ' a 



T + —dM^)) y e*™(2fc^x+^(x)) _ (2.1) 



kp 

n ■ 



Im Fall der Ladungsdichtewelle ist n = 2, da pro CDW-Periode zwei Elektronen vorhanden 
sind. 

In niedrigster Ordnung erhalt man folgenden Ausdruck 

p(x) = po + -dxfix) + pi cos(v?(x) + 2A;irx) + h.h. , (2.2) 
vr 

mit po = Ikp/TT und pi = 2po- 



2.2 Hamiltonian 

Die Hamiltonfunktion des Phasenfeldes ip setzt sich aus zwei Anteilen zusammen; 

Ti-kiass = T~iei + Himp, wobei Hei = f dx-^hvp {-g^'p)^ die elastische Energie und Hjmp die 

" Unordungsenergie" des Systems beschreibt (den Index j^jass lasse ich im weiteren Verlauf dieses 
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Klassisches Modell 



N 



Abschnitts weg). Das Unordungspotential istdurch V{x) = — ^ Vi6{x—Xi) gegeben, wobei Xi 

i=l 

die Position der i-ten Verunreinigung (/mpurity) und > dessen Starke ist. Dieses Potential 
koppelt an die Ladungsdichte p{x) und liefert drei zusatzliche Terme zum Hamiltonian: Einen 
konstanten Beitrag durch po, einen Gradiententerm und einen Beitrag durch die Kopplung an 
den schnell variierenden bzw. mikroskopischen Anteil der Ladungsdichte pi cos [(p{x) + Qx) 
mit dem Wellenvektor Q = Ikp der ungestorten Welle. Hohere Harmonische in der Dichte 
()2.2|1 werden nicht berucksichtigt. Der konstante Beitrag kann vernachlassigt werden und der 
Gradiententerm kann durch eine geeignete Umdefinition des Phasenfeldes in den elastischen 
Term absorbiert werden (siehe dazu auch im nachsten Abschnitt). 
Die Hamiltonfunktion des Systems der Lange L lautet daher |Gor77| : 



n 



dx 



hvp ( d 



Att V dx 



Vipi5{x - Xi) COS {(f + Qx) > 
i=i ) 



(2.3) 



Der mittlere Abstand der Verunreinigungen c-J^ = L/N ist gro(3 im Vergleich mit der CDW 
Wellenlange, also 





E 










^imp 


rji* 














Lc 


L 



Abbildung 5: Elastische Energie und mittlere Unordnungsenergie als Funktion der Lange 



Nun kann man eine Langenskala Lc definieren, die Fukuyama-Lee-Lange |FL78| . auf der die 
elastische Energie gleich der mittleren Unordnungsenergie wird {(p andert sich um O {!)): 
hvp/Lc = {V^piCimpLcY^'^ (siehe Abb. EI). Daraus folgt fur Lc 



Lc = {hVpfl\V^^p^C;^p)-^l\ 

Die durch Lc definierte Energieskala bezeichne ich als 

T* = ^^ = {nv,v!p.c-,^,f' 



(2.4) 



(2.5) 



2.2 Hamiltonian 
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Typische Werte fur T* /T liegen im Bereich 10^ - 10^ |IVIF91 |. 

Im Folgenden unterscheiden wir die Falle von schwacher {Vipi <C hvpcimp) und starker 
{ViPi » hvpCimp) Unordnung. 

Die angegebenen Bedingungen fur starke und schwache Unordnung kann man sich 
folgendermaBen Liberlegen: Im Fall eines starken Unordnungspotentials wird sich die Phase 
an den Positionen Xi der Impurities so einstellen, dass -ff/mp minimal wird. Dabei gewinnnt das 
System die Energie Vipi pro Verunreinigung. Andererseits wird dabei die elastische Energie 
erhoht, und zwar um ^ cJ^phvp^CimpY — ^^pCimp pro Verunreinigung. Die Unterscheidung 
von starker und schwacher Unordung ist dadurch vorgegeben, ob der Quotient Vipi/ {hvpCimp) 
sehr viel groBer bzw. kleiner als eins ist. 

Im Fall der schwachen Unordnung wird die Korrelationsfunktion {{'^{x) — (/?(0))^) 
zuerst mit einem einfachen Renormierungsgruppenargument und anschlieBend deren 
Temperaturabhangigkeit in one-/oop-Naherung mit einer KPZ-artigen Gleichung berechnet. 
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Klassisches Modell 



2.2.1 Schwaches Pinning - Burgers Gleichung 

Aus der Bedingung fur die schwache Unordnung und der Definition von Lc folgt, dass die 
Fukuyama-Lee-Lange groB im Vergleich zum mittleren Abstand der Impurities ist, also 

» c-^^ » . (2.6) 

Da wir uns fur die Physik auf Langenskalen > Lc interessieren, kann der Hamiltonian ()2.3|1 in 
Form eines XY-Modells mit Zufallsfeld geschrieben werden: 



7i 



wobei V"^ = V^piCjmp und a{x) eine Zufallsphase ist, mit a{xi) = IkpXi | 27r ( | steht fur 
modulo). FLir den Unordungsmittelwert von a{x) gilt 



e*("(^)-"(^'')) =5(x-x'). (2.8) 

In (|2.7j) wurde auBerdem beriicksichtigt, dass unter der Renormierungsgruppentransformation, 
welche von ()2.3j) auf ()2.7j) fLihrt, ein linearer Gradiententerm a / dx g{x)dx'p generiert wird 
und die elastische Konstante hvp nicht renormiert wird. 
Fur g{x) soil gelten 

g{x)=Q , g{x) g{x') = a5{x — x') . (2-9) 

Dass die elastische Konstante nicht renormiert wird, folgt aus der Tatsache, dass das Modell 
()2.3jl eine sogenannte statistische Tilt-Symmetrie |5VBQ88| besitzt. (In Anhang FA.S! ist dies 
ausfiihrlicher gezeigt.) 

Es ist zweckmaBig zu reskalierten GroBen uberzugehen: 

y = L^^x , L = L^^L 
(p{y) = Lp{x) , a{y) = a{x) 

9{y)=9{x) , a = Lea (2.10) 

rj-\ rj-\ l'^^ 

Damit laBt sich ()2.7jl umschreiben zu 

2 



'^=T-'ldy 



\ (If ~ ^^^0 ~ ^"'^ ~ ^^^^^ 



(2.11) 



AuBerdem brauchen wir das Hilfsfeld definiert durch 

d(j) _ dip 



ay -ay «(^' (^'i^' 



und damit 

y 



^{y) = m + J dz~g{z) + (^(0) - m . (2.13) 





2.2 Hamiltonian 
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In nullter Ordnung der Storungstheorie kann der nicht-lineare Term in ()2.11j) ignoriert werden 
und die Paarkorrelationsfunktion C{x) laBt sich in diesem Fall unter Benutzungvon (|2.13j) und 
()2.9|1 einfach berechnen: 



c{x) = {{^{x) - mf) = {{m - mf) 



4>iy)-m + / dzgiz) 





y 



(0(y)-0(O))") + J J dzdz' ~g{z)~g{z') 



= f|y|+^|y| = (^^ + (7) IxpC, (2.14) 
wobei fiir den Rauhigkeitsexponenten Q = 1/2 gilt. 

Die Berechnung der Korrelationsfunktion fiir das einfache elastische Modell (welche fur den 
vorletzten Schritt benutzt wurde) ist in Anliang rA.4l zu finden. 

Fiir T ^ hat Feigel'man [Fei^^l. ausgehend von ()2.3jl . mit einer Transfermatrix- 
Methode gezeigt, dass a von der GroBenordnung L^^ ist, d.h. a ist von der Ordnung 
1. Dies wird sehr gut durch die Monte-Carlo-Sinnulation in Abschnitt 12.2.31 bestatigt. Bin 
ahnliches Resultat wurde auch von Villain und Fernandez |VF84 | mit Hilfe einer Ortsraum- 
Renormierungsgruppenmethode gefunden. Mit ()2.14jl haben wir ein Resultat bekommen, 
welches auch bei endlichen Temperaturen gultig ist. 

Bei endlichen Temperaturen kann der EinfuB des nichtlinearen Terms der Unordnung mit 
einer KPZ- oder Surg'ers-ahnlichen Gleichnung behandelt werden. Dazu betrachtet man die 
eingeschrankte Zustandssumme 

Z{l,ip)^ I P^e-^«/^ (2.15) 
v(o)=o 

des Systems mit Lange / und den Randbedingungen 99(0) = 0, ip{l) = (p. 

Mit Hilfe der Transformation F{l,ip) = —Tln{Z{l,ip)) {Cole-Hopf-Transformation) erhalt 



man aus der "Transfermatix-Gleichung" _j'2^^j!£l = _|_ y(/^ (^)^ 2{l,ip) 

folgende Gleichung fur die eingeschrankte frei Energie F{l,ip) [HHF85j 



die 



Nun geht man zur Fouriertransformierten F{k,uj) = JJ^2 d(pdl e^^'^~^^'^'-F{l,(l)) (analog fiir 
V{k, u>)) Liber. Setzt man die Umkehrtransformierte in Gleichung ()2.1bj) ein, erhalt man [M+89| 
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T 

iujF{k,oj) = -— — k^F{k,uj) 
2hvp 



2 {k - k)kF{k -k,LO- Lo)F{k, Co) 



"^fiVF J J (27r) 

{[(i|<r2}xR 

+V{k,u). (2.17) 

Wird V{k,u}) explizit aufgeschrieben, sieht man, dass das Zufallspotential nur 
Fourierkomponenten mit = ±1 enthalt: 

V{k,uj) = V [[ d^dl cos{ip - a{l))e'^^+''^^ 
dip (^e<^+i)^ + e*(^-i)^ 



2 

h2{uj) 



+ 



dip 



2i 

= TiSik + l){hi{uj) - ih2{i^)) 

+T:5{k - l){hi{i^) + ih2{uj)) (2.18) 

dabei ist hi{uj) j cos(Q(x))e*'^' und /i2('^) analog mit sin(a) definiert. 

Damit kann jetzt das Unordnungs-gemittelte Produkt V{k' , iv')V{k, u) unter Ausnutzung von 

()2.8|1 berechnet werden: 



V{k',Lo')V{k,Lj) = 2TT^hi{Lo)hi{uj')6{k + k'){6{k - 1) + 6{k + 1)) 
= 2V^7r'^6{io + cj')'^(/fc + A;')(5(A; - 1) + 6{k + 1)) 
= 26{Lj + Lu')6{k + k')D{k,iv) (2.19) 

Nun fiihre ich folgende Bezeichnungen ein: 

(2.20) 



2hvp 



g(fe,cj)y(A;,a;) = F{k,uj) (2.22) 

' - ik 

ist im wesentlichen die inverse relativistische de Broglie-Wellenlange eines Teilchens mit 
Geschwindigkeit vp, Go der ungestorte und Q der effektive Propagator. 
Damit laBt sich Gl. (|2.17|1 als eine selbstkonsistente Gleichung fur Q{k,Lo) schreiben: 

g{k,iv)V{k,u;) = go{k,u;)V{k,iv)-^go{k,i^) II ^{k-k)k 

{|(i|<Q}xR 

xV{k -k,u- Co)g{k -k,iv- iJj)V{k,iv)g{k,Lv) (2.24) 
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G(k,o)) l/(k,co) Go(k,co) ^k,(o) Go(k,©; 




j l/(k-k,co-5) 
G(k-k,(o-S)) 



^k,co) 



Abbildung 6: Diagrammatische Darstellung der selbstkonsistenten Gleichung fur den effektiven 
Propagator (doppelte Linie). Der ungestorte Propagator ist durch eine einfache Linie und der 
Unordnungsterm durch ein x gekennzeichnet. 



Diagrammatisch ist diese Gleichung in Abb. Eldargestellt. Nun fuhrt man die Picard-lteration 
einmal durch, ersetzt Q durch Qq und erhalt so einen Korrekturterm der niedrigsten, nicht 
verschwindenden Ordnung (one-/oop-Naherung), welcher in Abb. [3a) dargestellt ist. Nach 
Mittelung uber die Unordnung, unter Benutzung von (|2.19j) . erhalt man 

GiKu:) = goik,u;)+i(^-^ygi{k,u;) jj - ~k) ■ k{-k) ■ k 

{\cj\<n}xM. 

xGoik -k,Lo- Cd)go{k,u)go{-k, -Cb)2D{k,Cj) + O(e^) , (2.25) 

siehe Abb.[7Ib). Der rechte Vertex entspricht dem Faktor {k — k) ■ k, der linke {—k) ■ k und 
O entspricht 2D{k,Lo). 

Fur u; = kann das Integral auf der rechten Seite von ()2.25j) weiter vereinfacht werden und 
aufgrund der einfachen Struktur von D{k, oj) ist das fc-lntegral einfach zu berechnen ^: 



g{k,o) = goik,ij) 



^ _ V^ire^k 



2 

i- J- 1 \ 1 1 

(2.26) 



k-1 k+1 

-T-, + 



v{k — 1)2 + lijj v{k + 1)2 + iCb J u"^ + oj"^ 

Um das a}-lntegral noch zu bestimmen, nutzt man aus, dass das physikalische Verhalten des 
Systems bei tiefen Temperaturen durch kleine k bestimmt wird und entwickelt den Integranden 
in (|2.26|) in eine Laurent-Reihe um A; = 0. Als relevanter Term bleibt nur noch ^^2^^"-)2 zu 
integrieren. FLir Q ^ oo erhalt man daher schlieBlich: 

g{k, 0) = go{k, u) (l - ^lv^k\^u-^\ (2.27) 



^Hat D{k,uj) eine kompliziertere fc-Abhanigkeit, z.B. D oc fc"^'' (raumlich korrelierte Unordnung), 
kann der Korrekturterm divergieren. In diesem Fall muss der Einfluss der Unordnung schrittweise mit einer 
Renormierungsgruppen-Rechnung behandelt werden ' M^89| . Fiihrt man diese Rechnung im vorliegenden Fall 
durch, erhalt man letzendlich dasselbe Ergebnis 112.2711 . 



14 



Klassisches Modell 




kjO) k-k,co-c5 k,co 

(b) 

Abbildung 7: Propagator in one-loop Naherung vor (a) und nach (b) Mittelung uber die 
Unordnung, wobei Q)=2D{k,u}). 

Vergleicht man die rechte Seite von (|T77ll mit g{k, 0) = ^ + findet man, dass die 

Temperatur oder besser die DifFusionskonstante Ihvp/T aufgrund der Unordnung in niedrigster 
Ordnung durch 

?^(l_4.= 5!|i^). (2.28) 
ersetzt wird. Der Rauhigkeitsexponent bleibt unverandert C, = 1/2. 

Die durch Gl. (|2.28|1 gegebene Abweichung der Temperaturabhanigkeit der Korrelationsfunktion 
vom ungestorten Fall, erfaBt aufgrund der one-loop- Naherung nur das Verhalten fur hohere 
Temperaturen T > (47r^)^/^T* bzw. kleine Korrekturen zur DifFusionskonstante. Daher 
liefert diese Naherung fiir groBe Temperaturen wieder das einfache lineare Verhalten, aber 
auch eine Erhohung der efFektiven Temperatur, aufgrund der Unordnung. In Abb. [HI ist die 
Temperaturabhanigkeit der Korrelationsfunktionen dargestellt. 

Dass die Paarkorrelationsfunktion bei T = mit Unordnung bei endlichen Abstanden [y > 
0) nicht null ist, ist klar, da die Phase an den Positionen der Verunreinigungen auf einen 
Wert fixiert wird, und somit die perfekte Ladungsordnung zerstort wird und damit auch die 
langreichweitige Ordnung. 
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Abbildung 8: Temperaturabhanigkeit der Korrelationsfunktion im ungestorten Fall (punktiert) 
und in one-/oop-Naherung (durchgezogen). Gestrichelt ist das Verhalten fur tiefe Temperaturen 
im Unordungsfall extrapoliert dargestellt. 
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2.2.2 Starkes Pinning 



FLir Starke Unordnung {ViPi ^ hvpCj^p) geht man zweckmaBigerweise zu einem effektiven, 
diskreten Modell uber, indem man das Phasenfeld (p{x) jeweils zwischen den Orten der 
Verunreinigungen in der Zustandssumme / e~^'<'^''/^ ausintegriert (Lc ist jetzt nicht 
mehr groB gegenuber c^p) [FeiSOj; hier exemplarisch zwischen </?i = ip{xi) und 93^+1 
(e = {xi+i-Xi)/My 



'P='Pi+l 



Vip e 



hv 

'JnT 



lim 



lim 

M^oo 



k=l 



^ M-i _fiup ( Vi,k+i-Vi,k~i ' 

i T— r , 47rT ^ ^\ 2e 



.g 27rgTe 



(^f ■Ai-^j+¥'?,o+'^?,M-2{¥'i,0'Pi,2+'Pi,M-2¥'i,A/)) 



(2.29) 



Dabei wurden folgende Bezeichnungen eingefiihrt: 
• 'fi,k = ^{xi + ke) 



5^ 



A. 





• • 5 










(2 





-1 







o\ 





2 





-1 






-1 





2 












-1 ■. 







-1 











2 





Vo 







-1 





2 J 



Das M-dimensionale GauBintegral kann leicht berechnet werden (siehe ()A.2.8j) ) und man 
erhalt: 



Vip e 



Re 



hvp (Vi+l-Vj) 



(2.30) 



mit einer Konstanten R. 

Daher hat der effektive Hamiltonian die Form 

N 



Heff= 2^ <^ — Z Z ^iP^ COs{ipi + QXi 



i=l 



in Xi - Xj+i 



(2.31) 



Da wir das Verhalten bei tiefen Temperaturen untersuchen wollen, nehmen wir weiter an, 
dass Vi ^ T. ipi wird, um die Energie zu minimieren, nur Werte annehmen, fur die gilt: 



(fi + Qxi = 2iTni mit nj E Z. Nun setze ich rii^i—rii = hi + 



_Q_ 



27rc,v 



und — Xi 



wobei fur die Xi gilt: < xi < ^2 < . . . < x^+i < L {L war die Systemlange). {[x] = k ist die 
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GauBklammerfunktion, welche ein A; G Z liefert, so daB | [x] — k\ minimal wird, insbesondere 
< 1/2.) Der effektive Hamiltonian lautet damit 



^ hv, ^ - ^ 



7 



i=l 



•-•imp 



(2.32) 



wobei 7 



27rc,> 



SO daB I7I < i ist. 



AuBerdem kann man schreiben: 



n-l 



(27rn„ - Qxn) + ^('/'i - V'i) - (27rni - Qxi) 



^ri-l 



i=2 

27rc 



/mp 



27r 



, j=i 



27r 



(2.33) 



Im folgenden ist es notig, die Statistik geordneter System zu benutzen - hier genauer gesagt 
die Statistik der Abstande der Verunreinigungen + c^^. In Anhang FA.SI ist die Herleitung 
der Verteilungsfunktion der Abstande allgemein dargestellt (Grundlage ist das Buch von 
David |Dav7Qj ). Sind die Positionen x,; der Impurities gleichwahrscheinlich zwischen und L 
verteilt, erhalt man die Wahrscheinlichkeitsdichtefunktion (pdf) oder kurz Verteilungsfunktion 
der ei, p{ei): 

p{ei) = c/^p(l - j{ei + (2.34) 

mit — < < L — (siehe Def. von e,). Im thermodynamischen Limes kann diese pdf 

gut durch 



-imp —c- f 

— g '-imp'^i 



< €{ < 00 



-'imp 



genahert werden. 

Damit konnen die n-ten Momente fiir n > 1 berechnet werden: 



00 



c 
-e 

1 e 



c d"- 



c 9" 



X=c 



A=c 



Substituiert man x = A/c und benutzt ()A.2.1|1 . erhalt man 



(2.35) 



(2.36) 



n! V - 



(2.37) 



fiir n > 1, ei = und fur den Korrelator EiEj = -^5ij. 

^imp 
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Nutzt man aus, dass hi = und hihj oc 6ij, folgt fur die Korrelationsfunktion mit (|2.33|) : 



47r^( h 



2tt 



n . 



(2.38) 



Bei niedrigen Temperaturen T <^ hvpcimp, nimmt {hi — ^ — 7) den minimalen Wert an, was 



genau der Fall ist, wenn hi 



2-K 



^ +7 



also 



2lT 



+ 7 



27r 



+ 7 



l/c,> 



27r 



2-K 



+ 7 



(2.39) 



Da Q ^ c/mp konnen wir den kleinen Parameter a = ^^^^ einfuhren. Unter Beachtung dass 
[x + n] = [x] + n fur n G Ij und der Substitution x = Cimpf-i + 1 laBt sich vereinfachen zu 

2 



Cf= / dxe"^' 



2a 



(2.40) 



Dieses Integral laBt sich exakt berechnen (siehe Anhang IA.5|1 und man erhalt folgenden 
Ausdruck: 



1 



1 



2a^ 2Q!sinh(Q) 
12 720 ^ > 



(2.41) 



wobei (sinha) ^ im letzten Schritt fur kleine a entwickelt wurde (Reihe fiir csch ist in (|A. 2.2(1 
aufgeschrieben). 



Zu bemerken ist, dass man unter der Annahme, dass ^+7 



1/2 ■ 



im Intervall [-1/2,1/2] 



1 



bekommen 



gleichwahrscheinlich verteilt ist, das einfache Resultat Cf = f_y2 
wiirde. 

Bei hohen Temperaturen T ^ hvpcimp kann man vernachlassigen, dass hi eine diskrete 
Variable ist, und erhalt in diesem Fall ^(27r/ij — Qei)^^ = -^r-^j—- 
Eine sinnvolle Interpolationsformel dieser beiden Ergebnisse ist durch 



vr 



^imp 



7TT^C^mp^ T 



180 Q2 f^yp, 



(2.42) 



gegebeb. Vergleicht man dies mit ()2.14jl sieht man, dass auch hier der Rauhigkeitsexponent 
unverandert C, = 1/2 ist. 



2.2.3 numerische Untersuchungen 

Dm das Resultat (|2.14|) fiir schwache Unordung numerisch zu iiberprufen, wurde ein Monte- 
Carlo Algorithmus verwendet. Ausgangspunkt war die diskrete Version des Hamiltonians ()2.11jl 
mit g{y) = 0, d.h. nicht renormiert, mit welchem die Paarkorrelationsfunktion berechnet wurde. 
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Bemerkung zum Algorithmus: Fur die Phase wurden freie Randbedingungen benutzt und 
das Erreichen des thermischen Gleichgewichtes wurde uberpruft, indem die Daten von 
Programmdurchlaufen mit mindestens dreimal langerer Equilibierungszeit verglichen wurden. 
FLir den Mittelungsprozess wurde die erste Halfe der Monte-Carlo-Schritte nicht beriicksichtigt. 




Abbildung 9: Die Abstandsabhangigkeit 
der Korrelationsfunktion C{y) fur 
verschiedene SystemgroBen: N = 100,200 
und 400; T = 1. Die Daten sind uber 
50-200 Unordnungskonfigurationen 
gemittelt. 



Abbildung 10: Die 

Temperaturabhangigkeit der Funktion 
C{y)/y - auch uber 50-200 
Unordnungskonfigurationen gemittelt. Die 
Gerade zeigt die Temperaturabhangigkeit 
der Korrelationsfunktion ohne Unordnung. 



Abb. zeigt die Korrelationsfunktion C{y) als Funktion des Abstandes y bei fixierter 
Temperatur fur verschiedene SystemgroBen. Wie man sieht, hat man fur alle SystemgroBen eine 
lineare Abhangigkeit vom Abstand y und bestatigt somit den Rauhigkeitsexponenten C, = 1/2. 
In Abb. [TOl ist die (interessantere) Temperaturabhangigkeit der Paarkorrelationsfunktion 
abbgebildet. Die Berechnung wurde mit einem System der GroBe N = 200 durchgefuhrt. 
Fur y — s- sieht man, dass C{y)/y gegen einen Wert 0(1), genauer gegen a 0.65, geht, 
was wiederum das Ergebnis ()2.14|) plausibel erscheinen laBt. AuBerdem sieht man, dass fur 
T < 0.3 die Unordnung dominiert und die Korrelationsfunktion temperaturunabhangig wird, 
wo hingegen fiir T > 1.3 die Unordnung irrelevant wird (siehe auch Abschnitt I4.7|l . 
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2.3 Dynamik 

In diesem Abschnitt wird das CDW-System im Fall eines externen elektischen Feldes untersucht. 
insbesondere wird das Kriechverhalten bei endlichen Temperaturen betrachtet. 



A 




Abbildung 11: "dynamische Regionen" der kollektiven CDW-Geschwindigkeit bzw. des CDW- 
Stroms: Depinning Region bei T = 0, Kriech-Region fur T* > T > und FlieB-Regime fiir 

E » Et. 

Wie schon fruher gesehen, erzeugt die Unordnung auf Langenskalen > Lc Energiebarrieren der 
Hohe T* . Bei T = und ohne externes Feld bleibt das Phasenfeld an den Verunreinigungen 
gepinnt. Legt man nun ein elektrisches Feld an das System an, bleibt dieser Zustand auch bis 
zu einem kritischen Schwellenwert des elektrischen Feldes Et, welches durch die Barrierenhohe 
vorgegeben ist, (bei T = !) erhalten. Wird die Feldstarke weiter erhoht, konnen die Barrieren 
uberwunden werden, sodass das Phasenfeld eine endliche mittlere Geschwindigkeit hat, d.h. 
ein Strom flieBt {Depinning Ubergang, siehe Abschnitt I2.3.2|) . In Abb. ist dies vereinfacht 
am Beispiel eines Teilchens veranschaulicht. 

Bei endlichen Temperaturen, insbesondere im Bereich < T < T* , wird man fur > 
immer einen Strom messen, da die thermischen Fluktuation der Ladungsdichtewelle eriauben 
die Barrieren zu uberwinden [Kriech-Bewegung, siehe Abschnitt I2.3.3j) . In Abb. [^sind die 
Regionen der moglichen dynamische Prozesse dargestellt. 

Fiir weitere Untersuchungen wird der Beitrag des externen Feldes zum Hamiltonian benotigt: 
Das elektrische Feld koppelt an den langsam variierenden Teil der Ladungsdichte oc dx^f - Der 
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Abbildung 12: Energielandschaft bei T = fiir ein Teilchen bei verschiedenen Werten von E. 



entsprechende Hamiltonian hat die Form 



He 



J dx Exdxf{x) = const. + J dx E{p{x) 



(2.43) 



Im nachsten Abschnitt wird die Bewegungsgleichung hergeleitet, welche Ausgangspunkt fiir 
die numerischen Untersuchunungen ist. 



2.3.1 Bewegungsgleichung und Kriechgesetz 

Als Ansatz fur die Bewegungsgleichung benutzen wir eine Langevin-Gleichung, gegeben durch 



dip 6H , 



(2.44) 



wobei 7 ein kinetischer Koeffizient und r]{x,t) ein gauBverteilter thermischer Rauschterm mit 



(77) = und 
{r]{x,t)7]{x',t')) = 2T-f 6{x - x') 6{t - t') 

ist. 

In Gleichung ()2.44j) ist H der voile Hamiltonian inklusive TCext- 

Fiihrt man zusatzlich zu (|2.11|) noch folgende Reskalierung der Zeit t ein 

rp* 

T = l-j-t, 
Lie 



(2.45) 



(2.46) 



^Bisher war die Wahl des Hilbertraumes bzw. welche Werte tp annehmen konnte nicht festgelegt, d.h. man 
konnte ip als Winkelvariabel £ [0, 27r] oder (/? G R wahlen. In Anwesenheit eines elektrischen Feldes ist der 
Hamiltonian allerdings nicht mehr invariant unter der Transformation + 2n7r, daher muB p alle Werte 

annehmen konnen (siehe |SZ88p . 
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erhalt man folgende reskalierte (uberdampfte) Bewegungsgleichung 

dip{y,T) 



^'P{y, t) + sm{ip{y, t) - a{y)) + E + fi{y, r) 



(2.47) 



Diese Gleichung enthalt als einzige Parameter f = T/T* und E = E / E* , wobei E* = T*/Lc 
{E* ist von der Ordnung des T = Depinning-Schwellenfeldes Et)- 

Zur numerischen Untersuchung dieser Gleichung wird eine diskretisierte Version der 
Bewegungsgleichung (|2.47|) benotigt: 



At 



{(fi+i - 2(pi + (pi^i) 
+ sm{Lpi- di) + E + Ti{i,T) , i 



1...N 



(2.48) 



Der erste Term ist der eindimensionale Gitter-Laplaceoperator und di ist eine Zufallsphase, 
welche im Intervall [0,27r[ gleichwahrscheinlich verteilt ist. Der letzte Term fi{i,T), das 
thermische Rauschen, ist durch (|2.45jl definiert. 

Die mittlere kollektive Kriechgeschwindigkeit ist durch v{E,f) = (^L^^ J dy dr(p{y,T) 
gegeben iGru94j, bzw. in diskreter Version: 



v{E,f) 



1 Aip, 



N 



At 



(2.49) 



Zu bemerken ist, dass die Kriechgeschwindigkeit proportional zum CDW-Strom j^dw ist. 
Die Bewegungsgleichung wurde numerisch mit einem speziellen Runge-Kutta Algorithmus, 
welcher fur stochastische System geeignet ist |GH81| . integriert, wobei periodische 
Randbedingungen benutzt wurden. Der Hauptteil des benutzten C-Programms ist in Anhang 
lA.TI zu finden. 



T=0 



f\ \ ■^"--^^ berechnete 
i \ \ Kurve 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



angefittete 
Kurve 



Abbildung 13: Stromverlauf fur groBe E 
mit und ohne Unordnung bei T = 0. 



Abbildung 14: Differenz der Strome aus 
Abb. ini Eine Funktion oc E^^ ist 
angepasst worden (gestrichelte Linie). 



FLir Temperaturen T ^ 1 sind die Energiebarrieren nicht mehr relevant, d.h. die 
" Energielandschaft" ist flach, und die Ladungdichtewelle fuhrt gemaB ip 'yE eine gedampfte 
Bewegung aus. Dasselbe gilt natiiriich auch fiir ^ 1. In Abb. El ist der Verlauf der 
Geschwindigkeit mit dem angelegten Feld fiir hohe Felder mit und ohne Unordnung aufgetragen 
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(bei T = 0). In Abb.EJist die Differenz dieser beiden Kurven aufgetragen und fur groBe E 
mit einer Funktion oc E^^ angefittet, d.h. v nahert sich mit 1/ E dem Verlauf ohne Unordnung 
bzw. bei hohen Temperaturen, was von S. ScheidI und V.M. Vinokur gezeigt wurde | 5V98| . 
Im umgekehrten Fall T ^ T* andern sich die Energiebarrieren auf der Skala L > Lc urn 
Betrage von der GroBenordnung E. Wenn man von ()2.43j) ausgeht, sieht man, dass aufgrund 
des externen Feldes die Energiebarrieren (|2.5|1 einen zusatzlichen Term proportional zu 

EL{ip^f^ ^El(^^Y^ (2.50) 

erhalten. Hierbei wurde benutzt, dass cr € O {L^^) und der Rauhigkeitsexponent C, = 1/2 ist. 
Die Energiebarrieren E^ haben infolgedessen die Form: 

/ r \ 1/2 

E%{L) ^ cbT'' ^ ceE { — \ L, (2.51) 

wobei ± sich jeweils auf die Bewegung parallel und entgegengesetzt dem angelegten Feld 
bezieht und cb und ce Konstanten von der Ordnung 1 sind. Ohne externes Feld sind die 
Barrieren E^, wie schon gesehen, von der Ordnung T*. 

Da fur L ^ Lc die Barrieren (parallel zum Feld) maximal werden, kann mit Hilfe des Arrhenius 
Gesetzes folgender Ausdruck fur die Kriechgeschwindigkeit aufgeschrieben werden: 

2 V 

= 7£e--(^VT),i,h(^) , (2.52) 

d.h. v{E,T) ist im wesentlichen die Differenz zwischen den Wahrscheinlichkeiten fur eine 
Bewegung parallel bzw. antiparallel zum Feld. Der Vorfaktor in ()2.52jl wurde so gewahit, dass 
man fur T > (T* + \E\Lc) das lineare Verhalten v{E) ^ erhalt. 

Man hat in einer Dimension demnach ein Kim-Anderson Verhalten. Dieses analytische 
Kriechgesetz ist vollig verschieden von dem nicht-analytischen in d > 2 ^NatOOj: 

Ki?)ocexp(-|(§)') , /x = ^ (2.53) 

(T^ und E^ sind Parameter, die von der effektiven Barrierenhohe und der typischen Langenskala 
des Problems abhangen.) In d = 2 folgt die Kriechgeschwindigkeit einem Potenz-Gesetz [TS92J. 
In den nachsten zwei Abschnitten werden die Ergebnisse der numerischen Untersuchung 
behandelt. 



2.3.2 Depinning bei T = 

Zur Uberprijfung des Algorithmus wurde zunachst der Depinning-Ubergang bei T = 
untersucht. Fur den Schwellenwert des externen Feldes wurde ein Wert von Et ^ 0.22 
gefunden {E* ist tatsachlich von der GroBenordnung E^)- Der kritische Exponent ^, 
definiert durch v{E) oc Q{E — Et){E — Et)^ (im kritischen Bereich), wurde durch eine 
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Abbildung 15: Depinning Bereich. Die 
durchgezogene Linie wurde mit dem in 
der doppelt-logarithmischen Auftragung 
ermittelten Wert fur ^ erstellt. 



Abbildung 16: doppelt-logarithmische 
Auftragung des Depinning Bereiches mit 
angepasster Gerade. 



Geradenanpassung in der doppelt-logarithmischen Auftragung ermittelt und hat den Wert 
^ = 0.57 ± 0.07. Dieser Wert stimmt mit dem in anderen Arbeiten von u.a. Sibani und 
Littelwood [SL90J oder auch Myers und Sethna [MS93| ermittelten gut Liberein. In Abb. 
undlSlsind die ermittelten Kurven aufgetragen. 



2.3.3 Kriechen 

In diesem Abschnitt werden wir die CDW-Geschwindigkeit bei endlichen Temperaturen 
untersuchen und Gleichung ()2.52j) fur die Kriechgeschwindigkeit uberprufen. 
Fur die Simulation des Kriechverhaltens wurde eine Simulationszeit von IOOOtq und At = 
0.05ro benutzt. Dreimal langer Zeiten haben das Resultat im wesentlichen nicht verandert. 
Um ()2.52jl nachzuweisen und die Parameter ce und cb zu bestimmen, wurde zunachst T 
festgehalten und E variiert. In diesem Fall wurde eine SystemgroBe von N = 5000 benutzt und 
Liber 50 Unordnungskonfigurationen {cn^i = 1 . . . A^} gemittelt. Das Resultat war wiederum 
stabil gegenuber groBeren Systemenlangen. In Abb.Elsind die Resultate dargestellt. Die gute 
Anpassung an eine Gerade unterstiitzt, dass durch (|2.52jl die Abhangigkeit vom externen Feld 
richtig erfaBt wird. Um den Wert von ce zu erhalten, wurde ein iteratives Minimal-Quadrat- 
Verfahren verwendet und damit ce = 2.5 ± 0.2 ermittelt. Damit wird auch die erwartete 
GroBenordnung dieser Konstante als richtig bestatigt. Um den berechneten Stromverlauf mit 
Experimenten von S. Zaitsev-Zotov zu vergleichen, tragt man die Daten doppelt-logarithmisch 
auf fAbb.lT^ und vergleicht dies mit dem in Fig. 1 in [ZZ93J dargestellten gemessenen Verlauf 
(Die Experimente wurden mit TaSs durchgefuhrt). Wie man sieht, stimmt das Verhalten 
qualitativ gut uberein, insbesondere im Gultigkeitsbereich von ()2.52j) (T < 1, £^ < 1 
bzw. im Experiment T < 20i^, E < lOOV/cm). In Abb. [121 wurden auch Kurven fur 
den Ubergangsbereich T > 1 vom Kriechverhalten zum FlieBen eingezeichnet, welche im 
Experiment fur Temperaturen T > 20K zu finden sind. AuBerdem sieht man, wenn man den 
Stromverlauf mit Fig. 3 in | ZZ93| vergleicht (und beachtet, dass £' = 1 in etwa E = lOOV/cm 
entspricht), dass ()2.52j) das Kriechverhalten sehr viel besser erfaBt, insbesondere fiir kleine E, 
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als das dort vorgeschlagene oc e ^"/^-Verhalten fur sehr kleine E (siehe Abb. TTI] bzw. inset 
in Fig. 3 in ZZ93| ;) oder das oc e-'=/(^''^)-Verhalten fur grdBere E (siehe Abb. [201 bzw. Fig. 3 
in LZZ93]). 

Um cb zu bestimmen, wurde diesmal die Temperaturabhangigkeit von v{E,T) fiir fixiertes 
E untersucht. Da hier eine Reduktion der Fehlerbalken sehr schwer zu erreichen war, 
wurde diesmal eine SystemgroBe von = 2000 gewahit und iiber mindestens 500 
Unordnungskonfigurationen gemittelt. Die Resultate sind in Abb. zu sehen. Mit dem 
ermittelten Wert fiir ce, wurde die Funktion v{x) = ^ye^'^^^ smh.{Bx) mit x = l/T und 
B = ceE an die berechneten Werte angefittet und damit cb = 0.35 ± 0.10 gefunden. 
Insgesamt wird das durch ()2.52j) bestimmte Kriechverhalten gut durch die numerischen 
Simulationen gestutzt. 

An Abb.I^ist auch zu erkennen, dass der Strom j^-dw fui' tiefe Temperaturen bei fixiertem E 
saturiert. Dies stimmt mit dem in Experimenten gemessenen Verhalten | ZZ93I [7Z94| uberein 
(siehe z.B. Abb. 2 in [7793| V 
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20 40 60 80 100 

A/'(T) sinh (cB^/f) 



Abbildung 17: Anhangigkeit der Kriechgeschwindigkeit vom externen Feld bei fixierter 
Temperatur T. f{T) = ^j-e"^ und M{T) ist eine Skalierungs-Funktion, um das lineare 
Verhalten besser darzustellen, mit folgenden Werte: A/'(0.10) = 1.0, A/'(0.12) = 2.5, 
Ar(0.14) = 6.5 und A/'(0.16) = 9.0. Fiir ce wurde 2.5 benutzt (siehe Text). 
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Abbildung 18: Die Temperaturabhangigkeit der Kriechgeschwindigkeit fiir verschiedene jeweils 
fixierte Werte von E, welche an den Kurven stehen. Fiir eine tibersichtlichere Darstellung 
wurden die Daten skaliert oder verschoben. Fur E = 0.08 ist 'iv{x) dargestellt, fiir E = 0.10 
und 0.15 wurden die Daten um 6-10~^ bzw. 10~^ nach oben verschoben. Typische Fehlerbalken 
sind eingezeichnet. 
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Abbildung 20: Grafik analog zu Fig 3 in [ZZ93] 
(Erklarung siehe Text) 



Abbildung 21: Grafik analog zu Fig 3 
(inset) in | ZZ93 |. 
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3 Erweiterung des Modells 

In diesem Kapitel wird das einfache klassische Modell aus dem vorherigen Kapitel zu einem 
quantenmechanischen erweitert. Hinzu kommen ein kinetischer Term und ein phase-slip-Jerm. 
In den folgenden Abschnitten wird zunachst der vollstandige Hamiltonoperator aufgeschrieben, 
die Wirkung aus der Pfadintegraldarstellung berechnet und schlieBlich der " Mechanismus" der 
phase-slips diskutiert. Zusatzlich wird die Paarkorrelationsfunktion des 1 + 1-dimensionalen 
harmonischen Modells, d.h. ohne Unordnung und phase slips, berechnet, da diese in Kapitel|3 
wichtig sein wird. 

3.1 quantenmechanischer Hamiltonian 

Das quantenmechanisch erweiterte Modell hat die Form 

= Hkin + 'Hklass + 'Hdisl (3.1) 

wobei Hkin der kinetische Anteil ist, mit 

L 

/■ ,,2 -, 

7ikin= dx—P^, P(x),ip(x') =-ihd(x-x'). (3.2) 
J 2c L J 



P ist der zu ip konjugierte Impulsoperator und = p entspricht einer Massendichte. 
ii.disi ist der neu hinzugekommene phase-slip-Jerm, definiert durch 



y j dx cos j dy P{y) | 



'Hdisi = X I dx cos \ — I dyP{y) j , (3.3) 

wobei A proportional zur Wahrscheinlichkeit der Ausbildung einer Dislokation und q eine 
zunachst beliebige Zahl ist. In Abschnitt 13.31 wird dieser Term ausfuhrlicher motiviert und 
diskutiert. 

Der gesamte, erweiterte Hamiltonian lautet also: 



2 / X \ 

n = J ^P' + |(9,^)2-ycos(p(^(x)+Q(x)) + Acos I ^ JdyP{y) 



(3.4) 



p ist ebenfalls eine zunachst beliebige Konstante. Die elastische Konstante c wird aufgrund der 
Quantenfluktuationen i.a. verschieden von dem Wert hvp/{2iT) im klassischen Fall sein. (Da 
im folgenden nur der Fall der schwachen Unordnung benotigt wird, ist der Unordnungsterm 
analog zu ()2.7|) aufgeschrieben.) 

3.2 Pfadintegraldarstellung 

Nun gehen wir iiber zur Wirkung; dafiir fiihre ich folgendes "duales" Feld 9 ein, gegeben durch: 

dj = -^P (3.5) 
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(siehe auch Kapitel Isj- Die Kommutatorrelation lautet damit: 6{x),ip{x') = ttQ{x — x'). 
Damit vereinfacht sich Hdisi zu X f dx cos{q9). Da im folgenden der Unordnungsterm und der 
Dislokationsterm getrennt behandelt werden, wird die Wirkung auch getrennt fur beide Terme 
berechnet. Beide Terme zusammen lassen sich i.a. nicht als geschlossener Ausdruck in einem 
der beiden Felder aufschreiben. 

Beginnen wir mit dem bekannten Unordnungsterm (A = 0): Die Wirkung ist in diesem Fall 
durch 

VipVee'^'^''^^ = I Vipe~^'^^^ (3.6) 



definiert, d.h. im Funktionalintegral wird das ^-Feld ausintegriert. C^:, ist gegeben durch 



4 = 



dx 



27r2c 



I TT 



(3.7) 



Da ji (fur y = oder A = 0) in einer einfachen (normal geordneten) Form vorliegt und TYq im 
Impulsraum im wesentlichen in eine Summe von harmonischen Oszillatoren zerfallt, laBt sich 
diese Form von L^p (oder auch Zg) einfach mit Hilfe des Pfadintegrals herleiten (siehe z.B. 
IM88]). 

Um das ^-Integral einfach berechnen zu konnen, bringt man auf folgende quadratische 
Form 



-1 



dx 



2 C^TT^ 



+ 2{dx'^f - Vcos{pip{x) + a{x)) 



(3.8) 



und kann mit der Substitution dxO = dxO + ^dr^ das GauBintegral im Impulsraum einfach 
berechnen. 

Letzendlich erhalt man: 



hl3 

Sip = I dr I dx 



L2 { (dxVi^, r)f + v-\dMx, T)f}-V cos{p^{x, r) + a(x))\ (3.9) 



(/3 = ((feB)T)-i). 

Auch hier kann man (fur T > 0) zu reskalierten GroBen x = xX^^ und f = T{h(3)~^ 
ubergehen, wobei Ay = hfiv die relativistische de 6rag'//e-Wellenlange eines Teilchens mit 
der Geschwindighkeit v ist: 



df 



L/Xt 



dx 



(3.10) 



—(3VXt cos{p(f{x, f) + a.[x)) 
Fur den Dislokationsterm {V = 0) geht man analog vor. Die Wirkung ist in diesem Fall durch 

VipVOe-^'^''^" = I VOe'^^/^ (3.11) 



I 
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definiert, also wird jetzt, im Gegensatz zum Unordnungsterm das 99-Feld ausintegriert. Lq ist 
durch 



dx h 



-1 



dx 



I'K^c 2, vr 



gegeben. Die entsprechende Wirkung lautet 



So 



hf3 







dr / dx 







27r2c 



{v\d^9)^ + (dref) + XcosiqO) 



oder in reskalierten Einheiten: 

1 L/Xt 



h 



df 



dx 



hv 
2^c 



{diOf + {dfOf + (3\\t cos(g^) 



(3.12) 



(3.13) 



(3.14) 



D.h. der Dislokationsterm wird durch ein 1 + 1-dimensionales Sinus-Gordon-ModeW im Feld 9 
beschrieben. Die Bedeutung des Feldes 9 im Zusammenhang mit dem Tomonaga-Luttinger- 
Modell wird kurz in Kapitel [Hlerklart. 

Nachdem nun die Wirkungen hergeleitet sind, werde ich den zusatzlichen phase-slip-Term 
etwas besser motivieren und herleiten. 



3.3 Phase-slips 

Der Term Hdisi wurde von Pokrovskii [BPT82J fur ein Modell benutzt, mit dem ein 
Phasenubergang von einer inkommensurablen zu einer kommensurablen Phase in einem 
anisotopen System mit Dislokationen untersucht wurde. Dieser Term kann auf zwei 
verschiedene Weisen verstanden werden. Zunachst kann der Term in einer diskreten Version 
aufgeschrieben werden, also z.B. 




(3.15) 



Man sieht, dass Tidisi sine Superposition (und damit hermitesch) von Produkten von 
Translationsoperatoren e ~ * ist, d.h. bis zum Teilchen Nx an der Positon x wird die Phase 
if (!) urn zbgvr verschoben. Diese Verschiebung kann auch als Einfuhrung eines qir-Vortex bei 
X verstanden werden (Sac99| . Am Ort des Vortex muB das Phasenfeld eine Singularitat haben, 
d.h. das Phasenfeld undefiniert sein und damit der Ordnungsparameter verschwinden, was 
zur Folge hat, dass die Amplitude auf null gedriickt wird. Man beriicksichtigt somit indirekt 
Amplitudenfluktuationen. 

Wozu dient dieser " Mechanismus" ? Dazu kann man sich folgendes Bild Liberlegen: Man 
betrachtet ein System mit einer Verunreinigung, welche etwas um die Position Xi (gauBisch) 
ausgeschmiert ist, und einem an das CDW-System angelegten externen elektrischen Feld. Die 
Phase ist zunachst auf einen Wert tpi mit ipi\2TT = Qxi "gepinnt". Nun kann die oben 
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von oben: 




X / 








/ 








\ 







Abbildung 22: Mechanismus der phase-slips. Im Bild oben ist eine (gepinnte) 
Phasenfeldkonfiguration im Potential eines Pinning-Zentrum und einem externen Feld 
eingezeichnet. Links unten ist eine Aufsicht der Konfiguration mit Ausbildung eines ±2tt- 
Vortex-Paares gezeigt. Aufgrund des Energiegewinns kann das Vortex-Paar aus dem System 
gedangt werden (rechts unten). 



beschriebene Ausbildung einer Dislokation dazu fuhren, dass das Phasenfeld be! Xi lokal 
"entspannt" wird und es sich weiterbewegen kann. In Abb [221 ist dies etwas anschaulicher 
dargestellt. 

Der Term kann auch auf eine andere Weise, wie folgt plausibel gemacht werden. Man geht jetzt 
von einem 1+1-dimensionalen (XY-)Modell Sxy = K/2 J df J dx [{df(p)'^ + (0^0)^] aus und 
nimmt die hier nicht berucksitigten Vortex-Konfigurationen (in Raum und Zeit) mit <f^ dcp = 
27rqj hinzu [qj € Z ist die Windungszahl des Vortex und Cj ein beliebiger geschlossener Weg 
urn den Vortex herum). Folgt man der Rechnung von Kosterlitz ^ KT73j . findet man, dass diese 
Vortices einen zusatzlichen Energiebeitrag, welcher einer Coulomb-WW in zwei Dimensionen 
entspricht, liefert. Man erhalt somit ein 1 + 1-dimensionales Coublomb-Gas-Modell. Jose, 
KadanofF, Kirkpatrick und Nelson | JKKN77] haben gezeigt, dass sich dieses Coulomb-Gas- 
Modell (mit q = 2) fur tiefe Temperaturen auf ein Sinus-Gordon-Modell der Form (vgl. mit 
(4.5) in deren Arbeit) 

J dx j df^ [i^fHx, f)f + {d^Hx, f))^] + 2yo cos(27rif(^) (3.16) 
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abbilden laBt. Dies kann mit So verglichen werden, indem man 6 = —0 und f = vt setzt; 



df / dx 



1 

2^ 



A 

H cosiqi: IgO) 

9 



(3.17) 



mit g = hv/c, X' = X/c (siehe auch nachstes Kapitel). Identifiziert man K = und 
X/g = 2yo sieht man, dass beide Modelle dieselbe Form haben (mit q = 2). Der Parameter 
yo kontrolliert die Ausbildung der Vortices und laBt diese fiir kleine Werte unwahrscheinlicher 
werden. D.h. fur unsere phase-slips, wie schon eingangs erwahnt, dass diese erst fur groBe 
Werte von A wahrscheinlich sind. 



3.4 harmonisches Modell 



In diesem Abschnitt werde ich die Berechnung der Paarkorrelationsfunktion 

(((^(x, r) - r'))^)g bzw. ((0(x, r) — 6'(x', r'))^)^ fur das harmonische Modell 
(y = A = 0) 



So 



dr / dx 







(3.18) 



(entsprechend in 9) darstellen. Da die Rechnung ahnlich zu der in IA.4I verlauft. beschranke 
ich mich auf wenige Zwischenschritte. Mit (p{x,T) = ^.(^/j) Z^fc n ^^^^^^'^"^Vfc.n ifk,n = 
J dx /dre"*('''^+'^"^V(a;,T), k = 27rm/L und = 2'Kn/{h(3) sind die Matsubara- 
Frequenzen) erhalt man 



So = ".it 

k,n 



J{k, n) 



hp 



(3.19) 



Zusatzlich fuhrt man den Quellterm Sqlj] = — J dx J dr j{x,T)ip{x,T) ein, womit die 
Zustandssumme lautet: 



-So-SQlj] 



n / / dafc,„ci5fc,„e"^^^"'."+'''-"^+W^'"=-"°'='"+'"='"^'='"^ (3.20) 



k,n 

mit ipk,n = ak,n + ibk,n und jk,n = oik,n + if3k,n- Die Korrelationsfunktion {tp{^)if{x)) q laBt 
sich damit als Funktionalableitung 

TJ-^{x' -x,t' -t) 



j=0 



T 



L{hp) 



k,n 



schreiben und man erhalt 



9, ,2 ^ ^ 1 _ „lLL)n{T-T')+lk{x-x') 

mx,r)-y,{x',r')f)o = ^Y.- 



cL(5 



n.k 



{vkY + 



(3.21) 



(3.22) 
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oder analog fur 9 



L{h(5)h 



n,k 



(3.23) 



FLir weitere Zwecke formen wir ()3.23)) etwas weiter um: 



27r^c L 
L{h(3)h2^ 



dk 



-A 



y 1 



tk{x-x') ^ COs{uJn{T - t')) 



(3.24) 



mit (|A.2.3p und (|A.2.5p : 



vrc 



(n/3)n 



dk 



h(3 
2v\k\ 



.,,, cosh((f -|r-r-|)Hfc|) ' 
sinh(Mi;|A;|) 



A 



i;^ / A; 



34 



Renormierungsgruppe 



4 Renormierungsgruppe 
4.1 Uberblick 

Das Ziel dieses Kapitels ist die Herleitung der Flussgleichungen fur die Parameter des 
Modells ()3.4jl . Mit der Renormierungsgruppentransformation (kurz: RG) werden Freiheitsgrade 
sukzessive aus dem System entfernt, d.h. man bekommt effektive Parameter bzw. einen 
effektiven Hamiltonian (oder Wirkung) auf groBeren Langenskalen. Man laBt den Hamiltonian 
durch den Hilbertraum der Hamiltonians "flieBen", indem sukzessive eine Partialspur bzw. ein 
RG-Operator auf den Hamiltonian angewendet und anschlieBend das System wieder auf die 
urspriingliche GroBe reskaliert wird. Von besonderem Interesse sind Fixpunkte der RG, d.h. 
Fixpunkte des RG-Operators, an welchen das entscheidende Verhalten des Systems abgelesen 
werden kann. 

Da fur die Renormierung unseres Modells eine anisotrope Impulsschalen-RG benutzt wurde, 
werde ich die relevanten Schritte kurz beschreiben (siehe z.B. auch |Car96 | oder |NatOO| ): 




Abbildung 23: " Impulsschale" in 1+1 Dimensionen 
(i) S wird in S' mittels der Transformation 

s^s' = 7^[cS] 

renormiert, wobei TZ wie folgt als Partialspur definiert ist 
mit der Zerlegung ip{x) = (p^{x) + (p^{x), wobei 

n 0<|fc|<A/fe 
A 

v>{x,t) = imLy^Y. E e<'^^+-"^Vfc,n (4.1) 

n \k\>A/b 



4.2 dimensionslose Parameter 
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(analog fur 9). D.h. die Partialspur wird uber Moden mit A/6 < |fc| < A gebildet; 
fur die Matsubara-Frequenzen ujn = 2TTn/{hf3) gibt es keinen Cutoff, d.h. dass das 
Feld If unabhangig von tj„ in der Impulsschale ausintegriert wird. 6 = e' ist der 
Reskalierungsparameter mit < / ^ 1. Der Cutoff fur \k\ ist durch A ?a a^^ gegeben. 
Die "Impulsschale" ist in Abb. veranschaulicht. 

(ii) Desweiteren mussen Reskalierungen von x, r und (/? vorgenommen werden, um die Dichte 
der Freiheitsgrade und der Fluktuationsstarke zu erhalten. Diese sind gegeben durch 

X ^ x' = x/b entsprechend k ^ k' = kb 

r T' = T/b^ uj^J = ub^ 

if ip'{x/b) = ipb-^ . (4.2) 



(iii) Um die Flussgleichungen der Parameter in S zu erhalten, nehmen wir an, dass der 
Parameter unter der RG (in erster Ordnung) zu 7^- = — /({7j})ln6) ubergeht, 
wobei /({7i}) eine Funktion von alien auftretenden Parametern sein kann. Berucksichtigt 
man noch die Reskalierung von jj: 7^-^ = 7^6" ^ 7j(l + alnb) erhalt man mit 7^- = 
7j + d'-fj folgende Flussgleichung: 

^ = (a-/({7.}))T,. (4.3) 



Bevor die RG fur unser Modell durchgefuhrt wird, werde ich die das System charakterisierenden 
dimensionslosen Parameter im nachsten Abschnitt einfuhren. Danach werden die Wirkungen 
Su, und Sg renormiert. 



4.2 dimensionslose Parameter 

Da die relativistische de Broglie-\Ne\\en\'ange At eine wichtige Langenskala ist, auf der 
Quantenfluktuationen relevant werden, kann man diese Lange mit der kleinsten auftretenden 
Lange A^^ ^ a vergleichen: 

Der so definierte Parameter i ist die reskalierte Temperatur und g kontrolliert 
die Quantenfluktuationen, da dieser Parameter im wesentlichen den Koeffizienten des 
Impulsoperators in p.4jl darstellt. 

Die dimensionslosen Parameter, welche die Unordnungsstarke und die Starke der Dislokationen 
beschreiben, sind im wesentlichen das Verhaltnis der nicht reskalierten Parameter zur 
elastischen Konstante. 
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Zusammengefasst lauten die relevanten Parameter unseres Modells 



t = —- (4.4) 

~g = h-^K'^-^ (4.5) 

V = (^A"'") ("-6) 
4.3 Replica-Wirkung / Unordnungsmittel 

Urn den Unordnungsterm zu renormieren, geht man im allgemeinen zum Replica-Hamiltonian 
bzw. hier zur Replica-Wirkung uber, d.h., dass der Unordnungsmittelwert ausgefiihrt wird. 
Der Replica-" Trick" besteht in folgender Darstellung der freien Energie 



1 



T = -TlnZ = -Tlim -{Z^ -1) . (4.8) 
n^O n 

und der Vertauschung der zwei Grenzwertprozesse (Unordnungsmittelung und Limes n — > 0)^, 
um die Mittelung des Logarithmus zu vermeiden: 

^^=11111^(^-1). (4.9) 
Die Replica-Wirkung 5"[{99'>'}^=i..n] wird definiert durch 



7=1-^ 




g-5"[{ip^}7 = l..n] _ (4-10) 



5" ist also eine Wirkung in der jeder Freiheitsgrad n-fach repliziert ist (a, 7 und v sind die 
Replica-Indices). Bei einigen Systemen kann die Vertauschung von Unordnungsmittelung und 
n ^ zu Problemen fuhren, in unserem Fall wird dies allerdings nicht angenommen. 
Fur die Wirkung ()3.9j) laBt sich die Replica-Version mit e^fc) = e^^^-^/^ fur gaul3verteiltes r 
leicht berechnen: 



g-5" ^ g-Ea'5o[</'"]g5E<,,.'5;^,p[¥''^,{a;,T}]5,„,p[¥>^{x^T'}] 



^ ^-Y.o.iSo[v"]-^T.,S dxdx' /dra!r'cos(vP"+a(x))cos(</,'"+Q(a;'))) ^4 

wobei <So der gauBsche Anteil von Sip (j3.18jl ist und Simp der Pinning-Anteil. 
Mit ()2.8j) kann der Unordnungsmittelwert leicht berechnet werden 



cos((/?"(x, r) + a{x)) cos((/?^(x', r') + a(j;')) = cos(v9°(x, r) - 99'' (x', t'))6{x - x') . (4.12) 



*which is "believed to create no trouble" 
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Die Replica-Wirkung zu p.9j) lautet also 

r-nP fL 



1 ^0 ^0 

aM=l 



dx 



((a,vp"(x,r))2 + ^;-2(9,vp-(x,r))2) 



(4.13) 



4.4 RG des Unordnungsterms 



In diesem Abschnitt werde ich die RG der Wirkung S^p bzw. der Replica-Version ()4.13j) etwas 
detailierter, insbesondere da dies bei endlichen Temperaturen geschieht, und allgemeiner in 
d+ 1-Dimensionen ausfuhren. 
Zunachst wird das Funktional 

R[f{D]^^cos{pf{r)) (4.14) 



definiert, womit der Unordnungsanteil der Replica-Wirkung ()4.13jl als 

hl3 



3"^ = - JJ drdr' J d'^x^ R[ip''{x,T) - ip'' {x,t')] (4.15) 







geschrieben werden kann. 

AuBerdem wird folgende zeitlich nichtlokale Korrelationsfunktion benotigt (siehe ()3.21|1 ) 

(v.(.,r)^(^,rO> = ^Ep)2^- (4-16) 
Die renormierte Wirkung wird durch 

^-sr\{^V^\ = I Y\{Vipl) e-'5"[{^<.^>}l (4.17) 

definiert, wobei das Feld tp^ in einen langwelligen {<^\) und kurzwelligen {<pl.) Anteil, wie in 

()4.1|) angegeben, aufgespalten wurde. 

Mit 



S-[y^}] 5o"<[{^^}]+4">[{^^}]+5S[{¥'J}] (4.18) 

f.a/3 



S14WI}\ = Y^l ^^y^/'i£^[(VxV'|fer))2 + z;-2(9,c^|(x,r))2] (4.19) 



a=l ■ 

erhalt man so fLir 

S'^iWl}] = SS,<[Wl}] + ln2o"< + ((e-^SK^J}] _ i^^^^^ , (4.20) 

wobei (^(^0[{{p^}]^'^ = In (^0[{ip^}]^ der Kumulantenmittelwert des Funktionals O 

bzgl. 5o,> ist. Diese Kumulante wird in eine Storungsreihe entwickelt, d.h. die Unordnung wird 

deshal 



als kleine Storung betrachtet, und deshalb fur kleine e = ^ entwickelt: 



V n -E^^^'^ (4-21) 
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mit Ai = ( ( ^ J dr dr' J d'^x cos{p{(p°' — (p'^')) ) ) .Da wir die Flussgleichungen fiir 

\V'-=^ //o,> 
kleine Unordnungsstarken herleiten wollen, beschranken wir uns auf die niedrigste Ordnung der 

Storungsreihe. Die renormierte Wirkung ist in 1. Ordnung in e gegeben durch: 

S'^iWl}] = 5o"< - / drdr' I d'x {{R[^^ - v.n))o,> • (4-22) 

^ V ' 

-5i,< 

Urn fortzufahren, definiere ich als Abkiirzung: ^p"'^ = + p"^^ mit 

H^l^ = H^l{^.r)-p-^{x,T') (4.23) 
Damit wird R\pf''^\ um p'^^ in kleinen entwickelt: 

= R\pr\ + ^vpr\H^r + \^'vpr\{^rf + o{{prf) (4.24) 

Aus der Definition von R folgt sofort: = (-1)™^?^™^/] und nach einer 

Umbennenung der Replicaindices, kann man schreiben: 

= R{pr\ + '^^{^rx^i + R^'yprx^i^pi? - ^w>) + o^pr?) (4.25) 

Und fur ((i?b°'1))o,> folgt 

=0 

+0 ((((v^>''^)^»o,>) (4.26) 

(Zur Berechnung der Kumulante im quadratischen Teil wurde das Wick-Theorem benutzt.) 
Obiger Ausdruck bedeutet, dass auBer < noch Korrekturen zu V und, wie wir gleich 
sehen werden, auch zu g liefert. Im folgenden wird nur noch der interessante Anteil 



untersucht. Betrachten wir zunachst den Teilausdruck R!'\p°!^\ ((</'>)^)o>- ist R!'\p'^^ 



V- ^-^> 1 t;2 /■> d'^k 1 



i?"b<'1 V'>>0,> (4.27) 

Betrachten wir zuna 
-:p2i?[¥'<1 und 

^'^> ^0'> " 1^1^ +u;2 ^Tc\ (2^ ^ (t;A;)2 + a;2 (4.28) 

wobei das ^ andeuten soil, dass die Summe bzw. das Integral nur uber Impulse mit A/6 < 
|A;| < A lauft. Die Frequenzsumme Ial3t sich mit ()A.2.3|1 ausrechnen ^ 



c (27r)'^ 7 A; V 2 

X/b 

— i^dA^-1 coth(ATA/2) ln(6) . (4.29) 
_c 



'da wir bei endlichen Temperaturen sind, darf hier nicht zum Integral J du> ubergegangen werden! 



4.4 RG des Unordnungsterms 
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Wir erhalten folgende Korrektur des Unordnungsparameters e: 



e^e' = e{l- p^~gKd coi\i{g / {2t)) In b) 



(4.30) 



Nun kommen wir zum zweiten Term von ()4.27|1 . welcher eine Korrektur zum Parameter g 
liefert. Dazu wird eine Gradientenentwicklung von durchgefuhrt, wozu man zu neuen 

Koordinaten ("Relativ- und Schwerpunktskoordinaten in r-Rlchtung") 

t + t' = 2u , n G [0, 

r - r = 2w , w£ [-11(3/2, h(3/2] 

iibergeht und 99"'" um u\n w entwickelt: 



ct,a 



(4.31) 



Die Replicaindices sind hier gleich, da nach ()4.16j) ein Kronecker-Delta in der 
Korrelationsfunktion auftritt. Mit (pl3T|) wird B!' um entwickelt: 



(4.32) 



Bemerkung: Ware (^"" eine zeitlich lokale Funktion, wiirde R"[(p'^] = R"[0] = const, von 
den langwelligen Anteilen des Feldes unabhangig sein, d.h. man wLirde nur eine unbedeutende 
Konstante als Korrektur bekommen, die auch in ()4.32j) als Term 0. Ordnung auftaucht, welche 
ich im weiteren nicht mehr beachte.o 

Der voile Korrekturterm lautet (bei der Integral-Transformation zu u und w wurde die Jacobi- 
Determinante von 1/4 berucksichtigt): 

hf3 hf3/2 



du 



dw / d'^xR""[0] 



du 



■{x, u)w 



-hl3/2 



k,n 



(4.33) 



Die Frequenzsumme kann nun mit ()A.2.5|1 und anschlieBend das tt;-lntegral mit ()A.2.6|1 
berechnet werden, wobei insbesondere die Bedingung in ()A.2.6|1 beachtet werden mu(3. 
Zu bermerken ist, dass der Wert der Frequenzsumme mit w auf dem Integrationsintervall 
exponetiell abfallt und somit die Entwicklung (|4.31|1 bzw. (|4.32|1 rechtfertigt. Letzendlich erhalt 



man mit R""[0] = ^p'^ {u r) 
hp 



dT j d'^xidr^'^f 



( 



2 hv'^c 



dk- 



-1 



A;4 



1 



XTk/2 



A/b 



smh{XTk/2 



(4.34) 



Der (Massen-)Koeffizient des kinetischen Terms wird in erster Ordnung wie folgt durch den 
Unordnungsterm renormiert: 



2nv'^ 



2hv^ 



2hv'^ 



1 



g/(2t) 

sinh(5/(2t) 



In 6 



(4.35) 



Zu bemerken ist, dass es zur elastischen Konstante in dieser Ordnung keine RG-Korrekturen 
durch die Unordnung gibt. Diese wurde man erst mit A2 bekommen. Die entgultigen 
Flussgleichungen werden in Abschnitt 14. bl bestimmt. 
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4.5 RG der "phase-slips" 

In diesem Abschnitt warden die Renormierungskorrekturen zu g, i und A durch den phase-slip 
Term hergeleitet. Gleichzeitig wird auch die RG eines 1 + 1-dimensionalen Sinus-Gordon- 
Modells bei endlichen Temperaturen behandelt, welches auch durch Se beschrieben wird. In 
diesem Abschnitt wird der Index g weggelassen. 

Analog zu (|4.21j) wird auch in diesem Fall der Kumulantenmittelwert des phase-slip-Terms in 
eine Storreihe in X/h entwickelt: 

k=l 



wobei 



'5a[{%}] = ^ / ^^ / • (4-37) 

In erster Ordnung {Bi) erhalt man eine einfache Korrektur zu A und die zweite Ordnung {B2) 
liefert Korrekturen zu 5 und i. 
Beginnen wir zunachst mit 

Bi = (J^j dr j dx cos{q9)^^ , (4.38) 

und setzen wieder 6 = 9^ + 0^. Mit ()A.2.10|1 erhalt man sofort: 

((cos(g6')))o > = cos(g6'<) ((cos(g6'>)))o > + sm(g6'<) ((sm(g6'>)))o^> 

= cos{qe<)e'^^'^>^o,> ^ (4.39) 

und mit (fT^ : 

= WvCf = »"'(AtA/2) in (, . (4.40) 
Der Parameter A wird in erster Ordnung durch 



A 



A' = A (^1 - ^ coth(5/(2t)) In 6^ 



(4.41) 



renormiert. 

FLir die Betrage aus der zweiten Ordnung ist die Rechnung etwas aufwendiger, da hier die 
einfache Gradientenentwickung, wie sie im vorigen Abschnitt benutzt wurde, nur divergente 
Beitrage liefern wurde. Der Koeffizient B2 ist gegeben durch 

B2 = ^j drdr' j dxdx' ((cos(g0(x, r)) cos(ge(x', r'))))^ > . (4.42) 

Mit Hilfe des Wick'schen Theorems und der Entwicklung aller cos-Terme in 0> vereinfacht 
sich dieser Ausdruck zu 

{{cos{qe) cos(g6l')))o > = {cos{qe) cos{q9'))^^^ - {cos{q9))Q^^ (cos(g6l'))o^> 

« qHm{q9^)sm{qe'^){9>e'^)^^^ 

= ^{cos{qi9^ - e'^)) - cos{q{9^ + 9'^))) {9^9'^)^ ^ (4.43) 



4.5 RG der "phase-slips" 
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mit 6' = 9{x',t'). Der Term cos(g(0< + ^<)) wird im weiteren nicht mehr beriicksichtigt, 
da dieser im Fall x = x' und r = r', die zweite Harmonische zur " Grundschwingung" wird 
und daher unter der RG weniger relevante Korrekturen liefert als der erste Term in (|4.43|1 
| JKKN77llRHU80| . 

Die Korrelationsfunktion in (|4.43jl ist mit ()3.25jl wiederum einfach zu erhalten: 



, cosh (1 - 2|T - T'|/(ri/3))ili) 



Nun betrachten wir den Term cos{q{6^ — 9'^)). Wie schon erwahnt liefert die einfache, naive 
Entwicklung dieses Cosinus-Terms bis zur zweiten Ordnung in A0< =6^ — 9'^ divergente 
Korrekturterme (siehe z.B. Diskussion in |God92l Kap. 4] von P. Nozieres). 
Um diesen Term besser zu erfassen, mussen wir eine (einfache) Operatorproduktentwicklung 
durchfuhren, welche alle Ordnungen der Cosinus-Reihe mit berucksichtigt. Dazu wird aus 
jedem Term der Form (A0<)^" fur n > 1 der relevante marginale Anteil (A6'<)^ extrahiert. 
Dies erreicht man, indem man zunachst folgenden Mittelwert betrachtet: 

(cos(gAe<))o^< = 1 + ({qA9^fp^ ^-^{qA9^f-~''^ (4.45) 
Jeden Summanden in dieser Darstellung kann man nun teilweise auswerten: 

Der zu untersuchende Term kann mit dieser Uberlegung geschrieben werden als 

cos(,A0<) . 1 + igA9,f Hj! ((,A^<)-->„^^^ , (4.47) 

d.h. fur Terme hoherer Ordnung wird jeweils ein Teilmittelwert ausgefuhrt, um den marginalen 
Anteil zu erhalten. Mit Hilfe des Wick-Theorems kann man 

ausrechnen | Kad79l fkdQ80| und bekommt nun als Entwicklung des Cosinus-Terms: 

cos(qA0<) = 1 - ^(gA^<)2e'^^^^^<^'^o.< . (4.49) 



Da Korrekturen zu t und g berechnet werden sollen, wird, ahnlich dem Vorgehen im vorigen 
Abschnitt, {qA9^)^ in neue Koordinanten 

2s = X + x' , 2r = X — x' 
2u = T + t' , 2w = T — t' 
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umgeschrieben (Jacobideterminante (1/4)^) und in den " Relativ"-Koordinaten r und w 
entwickelt. Als relevante Terme bleiben iibrig 



(<,A«<)2 = ,2 



ds 



du 



(4.50) 



Insgesamt erhalt man nach einigen Umformungen und den Umbennungen u t und s — > x 
den Korrekturterm der RG zu 5o[^<]: 



/ \ \ 2 4 
TTC I A\ q 



,-hl3 pL fhp/i poo 

^1^1 , /, . , N I dr / dx I dw / dr 
2vh\hJ smh(ATA/2) Vo Jo Jo Jo 



r I + [ w- 



dx 



dr 



■((6»<(x+r,T+«))-0<(a;-r,r-M)))2)^ ^ 



cos(2Ar)cosh ((l - ArA/2^ In 5 



(4.51) 



Die Korrelationsfunktion wurde schon in ()3.25j) berechnet und lautet fur T > mit der 
Substitution k = AtA;/2 und f = 4r/AT, w = 4w/{h(3): 

{{9<{x + r,T + w) - e^{x - r,T - w)f)^^^ 

= ^ / coth(fc) - cos(fA:) ~ ^ ^ = ^Ug/{2t)-r,w) (4.52) 

9 J rC y siiiii(/c) y fi' 

Da die Integrale uber r und w nicht weiter berechnet werden konnen, definiere ich folgende 
Funktionen: 

roo /"I 



1 



64 



fx{y,z) = dr dwr e 



cos(zr) cosh((l — jt^Dz) 



JO 
oo /■! 



friy,z) = L dr I dww e 



JO 



cos(zr) cosh((l — 



Damit lassen sich die renormierte elastische Konstante und Masse in S in kompakter Form 
schreiben als: 



2^c 



(2^ 



V2^2c 



7r3 A2 



4 t4 sinh(g/(2t)) 



/,(5,5/(2t))ln6 



(4.53) 



und 



27r2c 



27r2c 



27r2( 



3 \2 



vr^ A 



4 t4 sinh(5/(2t)) 



/,(5,5/(2t))ln6 



(4.54) 



woraus Betrage zu den im nachsten Abschnitt dargestellten Flussgleichungen fur g und t folgen. 
Vorher schauen wir uns die Korrekturen bei T = an. Dazu gehen wir zu ()4.51j) zuruck und 
bilden in diesem Ausdruck den Limes T — > und finden: 



4.6 Reskalierung und Flussgleichungen 
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Mit den Substitutionen k = k/A, f = 2Ar, w = 2Avw und (|4.55|1 fiir k statt A, laBt sich die 
Korrelationsfunktion (|3.25j) wie folgt schreiben: 

((0<(x + r,T + w) - 6'<(x - r, r - w))'^)^^ ^ 

= 1 (l- cosCkf)e-'^^] (4.56) 

9 Jo k ^ ) 

= ^ (ln(tZ;2 + f2)/2 + ^ + j^(f^u))) mit (4.57) 
* 9 

Jf. + Jf X 

mit 7 = 0.577216..., der Eulerschen Konstanterf' . Die Funktion ^{f,w) ist so gewahit, dass 
sie fiir groBe f und w verschwindet, d.h. das asymptotische Verhalten der Korrelationsfunktion 
ist, wie in zwei (raumlichen) Dimensionen, logarithmisch mit |r|. 
Mit den Definitionen 

-1 l-OO l-OO 2 _ 

fx{9) = Tt: dr dw r^e-^^^^'^'^Kos{f)e-'^ (4.58) 
16 Jo Jo 

-I POO POO 2^ 

/O(^) = — dr dww^e-^^^^'^'^hos{r)e-'^ mit (4.59) 
16 Jo Jo 

f^{w,f) = {ln{w^ + f^)/2 + -f + n{f,w)) (4.60) 
erhalt man bei T = folgende RG-Korrekturen: 



27r2cy \2tt'^cJ \2tt'^cJ \ 2 g 



2tt'^cJ \2tt'^cJ \2tt'^cJ \ 2 g 

4.6 Reskalierung und Flussgleichungen 

Wir konnen nun zunachst einmal die Flussgleichungen fur A = in d + 1 Dimensionen 
aufschreiben. Zuvor mul3 aber der Beitrag der Reskalierung (Punkt (ii) in dem eingangs 
beschriebenen Schema) bestimmt werden. Die Reskalierung von S'^ liefert: 

c' = c6'^+--2+2C ^ T' = Tb' 
Fiir die dimensionslosen Parameter folgt: 





T' 




d 


~t 

9 




V' 


= v'ic'r 



\ -1/2 

^ y^2-a!/2-2C 



Zur Herleitung der Identitat * wurde Mathematica benutzt. 
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Die Flussgleichungen lauten damit und den in Abschnitt 14.41 herereleiteten RG Korrekturen (s. 
Punkt (iii) im RG-Schema): 



di 
dlnb 

dg 
dlnb 

dV 
dlnb 



[2-d- 2C] i 



1 - d - 2C - \KdV^ ( 1 



~9/m \ 

sinh{g/{2i))) 



2 - d/2 - 2C - Ya^d coth(5/(2t)) 



V 



(4.63) 
~g (4.64) 

(4.65) 



Mit lim^j^oo coth(2;) = 1 kann man diese Gleichungen fur T = aufschreiben: 



dg 
dlnb 

dV 
dlnb 



1 - d - 2C - ^KdV^ 

2 

2 - d/2 - 2C - ^~gKd 



V 



(4.66) 
(4.67) 



Durch die Reskalierung von So erhalt man folgende Beziehungen {9' = Ob'^'): 



r 

2^ ' 



1 \ ^l-z+2C 



A' = A6i+^ 



Aus der Kommutatorrelation fur 9 und y?, welche auf alien Langenskalen die gleiche Form 
haben soil, folgt sofort C,' = —C, und damit fur die dimensionslosen Parameter folgt: 

t' = tV-2C , ~g'=~gb-^i , A' = A62-2C, (4.68) 

d.h. in d = 1 fiir t und g dieselben wie fur 5^. Mit den im vorigen Abschnitt hergeleiteten 
RG-Korrekturen erhalt man folgende Flussgleichungen fur i, g und A: 



di 
dlnb 

dg 
dlnb 

d~X 
dlnb 



1 - 2C 



7r3 A2 



4 t4 smh{g/{2t)) 



fr{g,g/i2i)) i 



(4.69) 



-2C 



7r3 A^ 

8 t4 sinh(5/(2t)) 



iUg,g/i2i)) + frig,g/{2i))) g (4.70) 



2 

2-2C-^coth(5/(2t))) A 



(4.71) 



Oder fijr T = 0: 



dg 
dlnb 

d\ 
dlnb 



-2C- 



9^ 



(4.72) 
(4.73) 



4.7 numerische Auswertung und Flussdiagramme 
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In den Flussgleichungen fur A wurden jeweils nur Korrekturen bis zur ersten Ordnung in A 
berucksichtigt. 



In d = 1 konnen die Flussgleichungen aller vier Parameter fur T > kombiniert werden und 
man erhalt zusammenfassend folgenden Satz von partiellen DifFerentialgleichungen: 



di 
din 6 

dg 
din 6 

dV 
din 6 

d~X 
din 6 



[1 - 2C - Fr] i 



aim \ 

smh{g/i2i))J 



{F, + Fr)/2 



^-2C-^gcotH~g/{2i)) 



V 



2-2C- ^g-' coth{g/{2i)) 



(4.74) 
~g (4.75) 

(4.76) 

(4.77) 



wobei Fx und F^- wie folgt definiert wurden 



F./r 



fx{y,z) 

fc{z;r,w) 
h{z; r, w) 



TT 



3^2 



4 t4 sinh(g/(2t)) 



fx/rig, g/C^i)) mit 



dr 

64 Jo JO 

^ />oo 

— / dr 

64 Jo JO 

dt 



tsinh(t) 

, . cosli((l — |w|)z) 
cos(zr)- 



dwr^e '^v ■^''^^'^'^^h{z;r,w) 
d'ww'^e~^^''^^'^'^^h{z]r,w) mit 
cosh(t) — cos(tr) cosh((l — |w|)t)] 



(4.78) 
(4.79) 
(4.80) 
(4.81) 



sinh(z) 

Korrekturen zu V bzgl. A und umgekehrt wurden nicht berucksichtigt, da dies Korrekturen 
min. 2. Ordnung waren und man daher davon ausgehen kann, dass diese nur kleine Beitrage 
liefern, andererseits sind diese Korrekturen im Hinblick auf den in Abschnitt l3.3l beschriebenen 
(moglichen) phase-slip-Mechanismus interessant. Allerdings ist aufgrund der unterschiedlichen 
Beschreibung der Unordnung und der phase-slips eine Berechnung (wahrscheinlich) sehr 
aufwendig. 

Zur Wahl von (: Die beiden Wirkungen 5^ und Se sind invariant unter dem Shift ip,6 —>■ 
{(f + 2mr), {6 + 2m:). Andererseits liefert die Reskalierung tp = (p'b^ im nichtlinearen Cosinus- 
Term einen zusatzlichen Term oc \n{b)C,<psm.{pip+a{x)) (analog fur 0, jeweils in erster Ordnung 
in In 6), sodass die reskalierten Wirkungen nicht mehr diese Symmetrie aufweisen. Daher sollte 
C = sein. 

4.7 numerische Auswertung und Flussdiagramme 

In diesem Abschnitt sollen die Flussgleichungen und Phasendiagramm ausfuhrlich untersucht 
werden. Dazu beginnen wir mit dem Fall A = 0. 
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Bisher war der Parameter p beliebig, urn moglichst allgemein zu bleiben. Um die 
Flussgleichungen aber weiter untersuchen zu konnen, mussp fixiert warden. Im Fall der CDW 
ist p = 1, da das Unordnungspotential an die Ladungsdichte p{x) koppelt und \n p p = 1 ist. 
AuBerdem wird C = gesetzt. 

Betrachten wird zunachst den einfachsten Fall T = in d = 1. 




Abbildung 24: Fluss von ^ fur verschiedene Abbildung 25: Fluss von V in der 

Startwerte in der delokalisierten Phase. delokalisierten Phase. 




Abbildung 26: Fluss von g in der Abbildung 27: Fluss von V in der 

lokalisierten Phase. lokalisierten Phase. 



Die Flussgleichungen werden mit einem Runge-Kutta-Verfahren numerisch integriert. Typische 
Fliisse der Parameter g und V sind in den Abbildungen I^^Tj?71 aufgezeichnet. 

Zeichnet man die Flusse in der g-V-Ebene erhalt man das in Abb. [2H1 numerisch bestimmte 
Phasendiagramm. Wie man sieht, gibt es zwei unterschiedliche Bereiche (Phasen) der 
Startwerte: eine lokalisierte Phase, in der die Unordnung relevant wird und V fur jeden 
Startwert aus der Phase gegen unendlich flieBt (und g gegen 0), und eine delokalisierte Phase, 
in der die Unordnung irrelevant wird. Beide Phasen werden durch eine Separatrix voneinander 
getrennt, welche im Fixpunkt = Svr, V = 0) endet. 

Diese einfachen Flusse in der g-V-Ehene, d.h. die Funktion V{g) kann auch analytisch 
bestimmt werden. Dazu werden die beiden partiellen Differentialgleichungen (PDGLn) 



4.7 numerische Auswertung und Flussdiagramme 
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Abbildung 28: Typische Flusse bei T = in der y-^-Ebene. Anfangswerte sind mit einem Kreis 
gekennzeichnet, die Flussrichtung mit einem Pfeil. Kreuze stehen fur Fixpunkte der Parameter. 
Die Separatix ist fett gezeichnet. 



durcheinander dividiert, womit man die einfache gewohnliche DGL 



dg V \ g J 
erhalt. Diese DGL hat fiir y > folgende Losung ^ 



(4^82) 



V(g) = ^2{g-3r[n{g)+c,), (4.83) 

wobei Cg eine Konstante ist, welche von den Anfangswerten abhangt. 
Fiir die Separatrix mit y(37r) = erhalt man folgende Gleichung: 

Vsep{g) = [2g-6n{l+ln{g/{3n)))]'/^ , g > 3n . (4.84) 

Entwickelt man diesen Ausdruck um den Fixpunkt, erhalt man einen kritischen Exponenten 
der Separatrix an diesem Fixpunkt von k = 1: 

Vsepig) {g - 37rr , n = l, (4.85) 

was auch gut in Abb. [2H1 (gestrichelte Linie) zu erkennen ist. 
'^einfach mit Separationsansatz integrierbar. 
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Vt 




0,0 371/2 471 671 871 f 

Abbildung 29: Typische Flusse in der i-V-Ebene. 

Nun gehen wir zu T > uber und betrachten zunachst den Fluss in der t-y-Ebene mit g = 0: 
In diesem Fall haben die DGLn wieder eine einfache Form und konnen analytisch gelost werden: 

= cie' (4.86) 

V{1) = c2e-"^^'/"+3/2' (4.87) 

oder 

V{i) = Q?/2g-i"M . (4.88) 

Der Wert der Konstanten q hangt wieder von den Anfangswerten ab. Typische Flusse sind 
in Abb. 1211 dargestellt. Man sieht, dass bei beliebig kleiner, aber endlicher Temperatur die 
Unordnung irrelevant wird. 

Dass dies auch bei endlichem g der Fall ist, wird in Abb. |3U1 und EDverdeutlicht, in denen der 
Fluss einmal in Ebenen mit konstantem g > projiziert ist, und in Abb. ED in Ebenen mit 
konstantem i. 

In Abb.|211erkennt man auBerdem, dass es eine ausgezeignete Temperatur ic = 37r/2 bei ^ = 
gibt, bei der sich der Fluss von V mit i umkehrt. 

FLir endliche g wird ic{g) (bzw. (*)) durch die Gleichung 3/2 — g/(27r) cot]i{g / {2ic{g))) = 
definiert. Die Losungskurve ic{g) ist in Abb. dargestellt. Im Bereich des Fixpunktes 
{Qci = 37r,t = 0) kann der Verlauf durch gci(t) oc (37r — i)"' mit dem kritischen Exponenten 
a = 0.250 ± 0.002 angenahert werden. 

In Abb. EH ist der Fluss aller drei Parameter (also $(/) = {i{l),g{l),V{l))) dreidimensional 
dargestellt. 



Bevor wir die Flussgleichungen fur F = 0, A > untersuchen konnen, muB der Verlauf der 
Funktionen fx ()4.79j) und ()4.80|) bestimmt werden. Dieserwurde numerisch fur verschiedene 
Werte von i und g bestimmt. In den Abbildungen B^TOTl ist der Graph dieser Funktionen 
dargestellt. Zwischen zwei benachbarten Werten, wurden die Funktionen kubisch interpoliert. 
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Abbildung 30: Flusse in t-F-Ebenen 
fiir verschiedene Werte von g > 0. 
An den Graphen stehen die Werte von 
g{0). Startpunkte sind mit einem Kreis 
gekennzeichnet. (Da hier der Fluss zu 
kleinen g nicht dargestellt wird, sieht es 
so aus, als ob V{i) immer bei i ^ 37r/2 
ein Maximum hat, was allerdings nicht der 
Fall ist - siehe Text.) 
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Abbildung 31: Flusse in ^-F-Ebenen fiir 
verschiedene Werte von t > 0. An den 
Graphen stehen die Werte von i{0). 



Mit diesen Funktionen werden im folgenden die Flussgleichungen numerisch untersucht. Da die 
numerische Berechnung dieser Funktionen allerdings sehr zeitaufwendig ist (insbesondere fiir 
groBe g und g/{2i)), ist hier eine numerische Untersuchung der Flussgleichungen nur begrenzt 
moglich. 

Bemerkung 1: Da die Funktion fc bei endlichen Temperaturen fur groBe f linear mit f ansteigt 
- d.h., dass die Paarkorrelationsfunktion sich aufgrund der beschrankten r-Ausdehnung des 
Systems fur groBe f quasi eindimensional verhalt - konvergiert das f-lntegral in fx und fr und 
rechtfertigt wiederum die Entwicklung fiir kleine f. Die Konvergenz der tD-lntegrale ist in jedem 
Fall, durch die exponentiell abfallenden Funktionen, garantiert. Bei T = ist die Konvergenz 
des f-lntegrals in (|4.58|1 wegen des logarithmischen Verhaltens der Paarkorrelationsfunktion 
(s. ()4.b0j) ) fiir bestimmte g-Werte nicht mehr gegeben, da das Integral 

oo 

Ia{a) = J f~"cos(f)(if , a>0 

a 

fiir a < nicht mehr existiert (mit linia^Q I a{a) = — sm(a)). Diese Bedingung definiert mit 
()4.58j) einen kritischen Wert fiir g, denn fijr groBe f gilt a « 2 — q'^n/A, womit nur fiir 
g < existiert. o 

Bemerkung 2: Wie schon in Abschnitt 13.31 und Abb. [221 angedeutet, setzen wir q = 2. 
Diese Wahl ist plausibel, da zum einen die Werte von q € 2Z\{0} sein miissen (siehe Abb. 
r22|) und auBerdem q den minimalen Wert annimmt, da die Wahrscheinlichkeit fiir phase-slips 
if —>^ if + qn fiir groBere g = 4, 6, . . . bei tiefen Temperaturen sehr klein sein wird. o 
Kommen wir jetzt zur Untersuchung der Flusse bei T = 0. An den Flussgleichungen fiir A 
()4.77j) sieht man, dass es wieder eine Linie gcii^) gibt, bei der sich diesmal der Fluss von 



50 



Renormierungsgruppe 




g 



Abbildung 32: kritische Linie in der V = Ebene. Gestrichelt ist die Kurve ic{g) oc (Svr — 
mit a = 0.250 ± 0.02 am Fixpunkt {g = Svr,^ = 0) angefittet. Die Pfeile sollen den trivialen 
Fluss in dieser Ebene verdeutlichen. 



A umkehrt. Diese ist durch die Gleichung 2 — tt / gc2{i) coth.{gc2{t) / {2i)) = definiert. Fur 
i = erhalt man einen kritischen Wert von gc2 = 7r/2. Der Verlauf ist in Abb. BHldargestellt. 
Diese Werte von g trennen also, ahnlich wie die Kurve ic{g) fur V, zwei Regionen in denen 
A gegen null {g < 5c2(^)) bzw. gegen unendlich [g > gc2{i)) flieBt. Der Wert von stimmt 
mit dem Wert von g uberein, bei dem /]? (unbestimmt) divergiert, was bedeutet, dass die 
Flussgleichung fur g ()4.72|) bei T = nur bis g = 7r/2 gultig ist. Fur g > 7r/2 ist die gezeigte 
RG fur g nicht mehr ausreichend. Um dennoch Aussagen fur g > 7r/2 bei T = machen zu 
konnen, vernachlassigen wir in diesem Bereich die RG-Korrekturen bzgl. A zu g, was bedeutet, 
dass g{l) nicht mehr "flie[3t" {dg/dl = 0) und somit A fur alle g > tt/2 einfach gegen unendlich 
flieBt (sofern A(0) > war). 

Der komplette Fluss bei T = ist in Abb. EHl dargestellt. Man findet wiederum zwei Phasen 
(A und B), in denen einmal A gegen flieBt (A), d.h. die phase-slips irrelevant werden, und 
in der Phase B die phase-slips relevant werden. Da + zwischen und tt/2 bei 1.07 
eine Nullstelle besitzt (s. Abb. I34|l. ist die Steigung der Separatrix der beiden Phasen negativ. 
Aufgrund der Divergenz von f^{g) am Fixpunkt (5 = 5c2. = 0). ist ein kritischer Exponent 
der Separatrix am Fixpunkt nicht bestimmbar. Untersucht man allerdings die in Abb. 15^ fett 
dargestellte Flusslinie, welche naherungsweise die Separatrix ist, findet man, dass die Separatrix 
naherungsweise linear in den Fixpunkt verlauft. 

Der Fluss fiir endliche i ist in Abb. EHloder auch in Abb. |^ in Ebenen mit konstantem i 
projiziert, dargestellt. In Abb. |^ (rechts) und fiir die zwei Flusslinien in Abb. EHl bei groBen 
g ist zu beachten, dass die Differentialgleichungen, aus den anfanglich genannten Griinden, 
nur bis zu einem bestimmten Wert von g integriert werden konnten und nicht klar ist (aber 
wahrscheinlich, siehe nachstes Kapitel), ob der Fluss, wie die zwei "mittleren" Flusslinien in 
Abb. EOl die kritische Linie gc2{i) kreuzen und somit A letztendlich auch gegen null flieBt, 
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Abbildung 33: Dreidimensionale Darstellung des Flusses mit eingezeichneter Separatrix und 
icig) (gestrichelt). 

wie dies bei kleineren Startwerten von g{l) der Fall ist. Dies bedeutet also, dass auch die 
phase-slips, zumindest fur kleine Ausgangswerte von g, bei endlichen Temperaturen irrelevant 
werden. 



Zum SchluB schauen wir uns noch den kompletten Fluss bei T = fur V, A und g an, 
welcher in Abb. |421 visualisiert ist. Man kann drei Bereiche identifizieren: Fur g < gc^ ist 
die Unordnung relevant und die phase-slips irrelevant, im Zwischenbereich gc2 ^ 9 ^ 9ci 
sind beide Fluktuationen relevant (wenn A = oder V = 0) und den Bereich g > gc^ in 
dem die phase-slips relevant werden, aber die Unordnung irrelevant wird. Zur Berechnung der 
Flusslinien wurde fur g > gc2 fx — ^ gesetzt. 

Da in den Flussgleichungen fijr A und V keine RG-Korrekturen der jeweils andere Parameter 
berLicksichtigt wurden, sind die Aussagen fur die FILisse bei endlichem A und V mit Vorsicht zu 
geniessen. In den zwei Bereichen g < gc^ und g > gc^ sollte der Verlauf der Flusslinien durch 
die beiden "auBeren" Flusslinien in Abb. |421 richtig erfasst werden, da diese eine sinnvolle 
Erweiterung der bisher gefundenen Flusse sind. Im "Zwischenbereich" wird g zu kleineren 
Werten flieBen und V damit relevant bleiben, A kann dann aber wieder in den " irrelevanten 
Bereich" flieBen, was durch den Verlauf der berechneten Flusslinie unterstutzt wird ("mittlere" 
Flusslinie in Abb. I42j). Allerdings scheint dieser Fluss physikalisch nicht sinnvoll zu sein, was 
im nachsten Kapitel mit einem einfachen Argument dargestellt wird. 
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Abbildung 34: Funktionsverlauf von fx mit g 
bei verschiedenen Temperaturen i. Der Wert 
der Nullstelle dieser Funktion geht fiir t — > 
gegen den Wert ^ 0.92 (NS von /°) und wird 
fiir groBere i groBer. 




Abbildung 36: Funktionsverlauf von mit g 
bei verschiedenen Temperaturen i. fr{g = 0, t) 
wachst mit i. 




Abbildung 35: Funktionsverlauf von 

fx mit g und i. 




Abbildung 37: Funktionsverlauf von 

/t- mit g und i. 
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Abbildung 39: Fluss bei t = und V = 0. In der Phase A werden die phase-slips irrelevant 
und in B relevant. Die (naherungsweise) Separatrix ist fett dargestellt. 
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Abbildung 40: Flusse fur V = 0. Gestrichelt ist die Separatrix der beiden Phasen A und B und 
die kritische Linie gc2(^) eingezeichnet. 




Abbildung 41: Fluss fiir F = und T > in Ebenen mit konstantem i projiziert fur ^(0) < 
(links) und ^(0) > (rechts). 
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Abbildung 42: Fluss bei t = 0. 
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5 Resultate und Relevanz fiir verwandte Systeme 

In diesem Kapitel mochte ich die physikalischen Konsequenzen der Resultate aus dem vorigen 
Kapitel etwas ausfiihrlicher besprechen und teilweise auch mit denen schon aus hoheren 
Dimensionen oder anderen Modellen bekannten vergleichen. 

Bel r = im Fall ohne Dislokationen haben wir fur zunehmende Quantenfluktuationen 
einen Phasenubergang von einer lokalisierten, in der die Unordnung unter der RG wachst, 
zu einer delokalisierten Phase, in der die Unordnung zu null renormiert wird, gefunden. Dieser 
Phasenubergang findet bei einem Wert von = Svr/p^ in einer Dimension statt, wobei die 
Grenze der Phasen sowohl von der anfanglichen Unordnungsstarke, als auch von der Starke 
der Quantenfluktuationen abhangt (siehe Abb. I28|) . Dieses Phasendiagramm wurde schon 
von Giamarchi und Schuiz [GS88J fiir ein eindimensionales wechselwirkendes Elektronengas 
im Zufallspotential gefunden. In diesem Paper entspricht der Parameter K, welcher die 
Wechselwirkungsstarke kontrolliert (siehe auch Abschnitt 15. 2j) . dem Parameter g^^ bei uns 
(bis auf einen numerischen Vorfaktor). 

Nun konnen wir uns die Relevanz dieses Phasenubergangs fiir Ladungsdichtewellen anschauen. 
Der Hamiltonian ohne Unordnung hatte die Form 



Fiir CDWs ist c « ^ (siehe Kapitel \^ v ^ vf{^) ' |Gru94| . wobei nie die 
Elektronenmasse und m* die effektive Masse ist. FLir bekannte Ladungsdichtewellen-Systeme 
ist m* /nie von der GroBenordnung 10^ [GruBBj. Vergleicht man nun die GroBenordnungen 



dxff) haben dieselbe GroBenordnung) findet man, dass der kinetische Term um einen Faktor 
~ 10^ kleiner ist als der elastische Term, was die klassische Behandlung des Modellls in Kapitel 
121 rechtfertigt. FLir den Wert des Parameters g heiBt das daher, dass dieser in realisierbaren 
CDW-Systemen im Bereich von 10~^ — 10"^ liegt, also wesentlich kleiner als gc^ ist. Dies 
stimmt mit der bisherigen Beobachtung uberein, dass die Unordung in Ladungsdichtewellen 
immer relevant, d.h. die delokalisierte Phase nicht realisiert ist. 

Betrachtet man den phase-slip Term bzw. das entsprechende Phasendiagramm, findet man 
ebenfalls einen Phasenubergang bei T = bei einem Wert von = (Z^tt/S von einer Phase, 
in der die phase-slips irrelevant sind (A), zu einer Phase, in der diese relevant sind (B), siehe 
Abb. EHl Die Phasengrenze hangt wiederum von den Startwerten der Quantenfluktuationen 
und der phase-slip Wahrscheinlichkeiten ab. 

Aufgrund der Struktur von Sq, ist ein Vergleich mit dem Phasendiagramm eines 
zweidimensionalen Sinus-Gordon Modells bei T = angebracht (siehe z.B. [Sac99| ). Dieses 
Modell hat zwei Phasen, in denen der nichtlineare Term relevant bzw. irrelevant wird. In 
letzterem Fall flieBt der Koeffizient auf die Fixpunktlinie A = 0, 5 < ^C2. d.h. auf groBen 
Langenskalen kann das System durch eine freie GauBsche Theorie beschrieben werden. Das 
Phasendiagramm stimmt demnach (qualitativ) mit unserem uberein. Da, wie in Abschnitt 
13.31 angedeutet. ein XY-Modell mit Vortices auf ein Sinus-Gordon-Modell abgebildet werden 






unter Beachtung der Kommutatorrelation {P/h und 



5.1 Vergleich mit Flussen der Unordnung in d > 1 
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kann, laBt sich dieses Phasendiagramm auch mit einem Kosterlitz-Thouless Ubergang |Kos74| 
vergleichen. 

In unserem Fall hat dieser Ubergang entscheidene Konsequenzen auf den Lokalisierung- 
Delokalisierungs-Ubergang, da 

7r/2 = 5c2 < 5ci = 37r. 

Im Bereich g < wird in jedem Fall die Unordnung relevant und die Dislokationen 
irrelevant sein, und fiir g > gc^ dominieren die phase-slips. Allerdings ist das Verhalten 
des Systems auf groBen Langenskalen im Zwischenbereich gc2 < 9 < 9ci unklar. Die 
numerischen Ergebnisse sind in diesem Bereich, aufgrund der im vorigen Kapitel geschilderten 
Schwierigkeiten, nicht sehr verlasslich. Die Untersuchung lieferte, zumindest im Gebiet in 
dem die Flussgleichungen gelten, dass die Dislokationen in diesem Zwischenbereich irrelevant 
werden, was allerdings physikalisch nicht sehr plausibel ist, wie man sich folgendermaBen 
iiberlegen kann: Angenommen wir haben ein System in welchem phase-slips relevant sind, 
d.h. die Korrelationen in der Phase durch die Sprunge zerstort werden. "Schaltet" man 
nun schwache Unordnung ein, wird dies keinen Einfluss auf die PhasensprLinge haben, und 
diese auch auf groBen Langenskalen nicht unterdrucken. Daher ist es wahrscheinlich, dass der 
Lokalisierungs- Delokalisierungs-Ubergang bei verschwindet und der Ubergangswert von 
g um immerhin einen Faktor 6 zu gc2 vermindert wird. Fiir reale CDW-Systeme wird dies 
wiederum keine Konsequenzen haben, da wie schon erwahnt, realistische Werte fur g kleiner 
als gc2 sind. 

Bei endlichen Temperaturen werden sowohl die Unordnungsfluktuationen (siehe Abb. I33|l als 
auch die phase-slips irrelevant, sodass nur thermische Fluktuationen auf groBen Langenskalen 
eine Rolle spielen. Im Fall der Unordnung ist dies klar, da der coth-Term in der Flussgleichung 
fur groBe i linear mit i verlauft und somit, aufgrund des negativen Vorzeichens, V gegen null 
geht. Dieses Verhalten haben wir schon in Abschnitt l2.2l anhand der Temperaturabhangigkeit 
der Korrelationsfunktion gesehen. 

Fiir die phase-slips kann dies folgendermaBen erklart werden: da der Fluss von g fiir groBe i 
sehr stark abnimmt (siehe Def. von Fx, Ft und den Verlauf von fx und mit i), wird 
irgendwann auf jeden Fall die kritische Linie 5c2(^) (bei entsprechend groBen Startwerten von 
g) kreuzen und somit gegen null flieBen. Hierbei muss aber beriicksichtigt werden, dass die 
Flussgleichungen nur fiir kleine A giiltig sind. 

5.1 Vergleich mit Flussen der Unordnung \n d > 1 

Den Fluss in der i-V-Ebene ohne Dislokationen (A = 0) und ohne Quantenfluktuationen 
[g = 0), d.h. die Relevanz der Unordnung in eindimensionalen elastischen Medien bei endlichen 
Temperaturen, konnen wir mit dem Fluss in d > 2 vergleichen. 

In vier oder mehr Dimensionen werden sowohl Unordnungs- als auch thermische Fluktuationen 
irrelevant und man hat demnach einen Fixpunkt bei {T = 0, V = 0). Die 
Paarkorrelationsfunktion geht gegen einen konstanten Wert, d.h. es gibt keine Rauhigkeit. 
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Abbildung 43: Karikatur der Fliisse in d Dimensionen in der T-y-Ebene. 



In Dimensionen 2 < d < 4 bleiben thermische Fluktuationen irrelevant, der Fixpunkt "wandert" 
zu einem endlichen Wert von V (wobei zu beachten ist, dass hier funktionale RG benutzt 
werden muss), welcher in geeigneten Einheiten von der Ordnung ist. Die Rauhigkeit ist 
hier logarithmisch. 

In d = 2 gibt es eine endliche Ubergangstemperatur Tg (Glasubergangstemperatur) oberhalb 
welcher die Unordnung irrelevant wird, und unterhalb welcher sie auf eine Fixpunktlinie 
V* oc (Tg — T) flieBt. Die Glasphase ist durch ein superrauhes Verhalten gekennzeichnet 
(die Paarkorrelationsfunktion ist oc (Inr)^). Siehe dazu auch [NSOOJ . 

Die Flusse \n d > A, 2 < d < A, d = 2 und in einer Dimension, in welcher, wie im vorigen Kapitel 
gesehen, die thermischen Fluktuationen relevant und die Unordnungsfluktuationen irrelevant 
werden, sind in Abb. l4^zusammengestellt. 

5.2 Tomonaga-Luttinger Modell 

Aufgrund der Formulierung der Unordnung und der phase-slips in zwei verschiedenen Feldern, 
bietet sich ein Vergleich mit dem Tomonaga-Luttinger (TL) Modell an, mit welchem man (ohne 
Unordnung und Vortices) eindimensionale Elektronensysteme bei tiefen Temperaturen {ksT <C 



5.2 Tomonaga-Luttinger Modell 
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ep) beschreiben kann, in welchen die Dispersionsrelation an der Fermikante linearisiert werden 
kann. Der zunachst fermionische Hamiltonian kann in diesem Fall mit Hilfe der Bosonisierung 
(auf die ich hier nicht weiter eingehe) auf ein TL Modell abgebildet werden. Die Wirkung 
wird meist in einem von zwei zueinander "dualen" Feldern beschrieben, wobei das jeweils 
andere ausintegriert wurde - je nachdem fur welche Korrelationen man sich interessiert. Die 
autretenden Wirkungen sind auBerdem die gleichen, wie fiir ein masseloses, skales Feld. 
Identifiziert man g/n mit dem Luttingerparameter K, welcher ein MaB fiir die Dichte-Dichte- 
Wechselwirkung ist {K = 1: keine WW, K > 1: attraktive WW, K < I: repulsive WW), und 
V mit der renormalisierten Fermi-Geschwindigkeit, kann man S,pfi und Sq^ schreiben als 



d.h. als dieselbe Wirkung, wenn man K gegen 1/K austauscht. Das Feld ip hat im TL Modell 
die Bedeutung eines " phononischen" Verschiebungsfeldes und 6 ist das zu "duale" Feld, 
hat aber keine direkte physikalische Bedeutung. Der Gradient von 6, also der zu ip konjugierte 
Impuls, misst die Differenz in der Dichte von links- und rechts- laufenden Teilchen |5ac99| . 
also den Strom (vgl. dx^p oc p). 

Im TL Modell wird der Einfluss der Unordnung (insbes. backscattering) ebenfalls durch einen 
Cosinus-Term mit zufalliger Phase beschrieben (s. z.B. [GS88]). Vortices im Feld tp konnen 
auch durch einen Cosinus-Term im Feld 6 (oder umgekehrt) beschrieben werden | 5ac99 |. 




bzw. 
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6 Zusammenfassung und Ausblicke 

Ladungsdichtewellen-Systeme werden spatestens seit den 1970er Jahren intensiv theoretisch 
untersucht. Das auch heute noch neue interessante Resultate und Ergebnisse gefunden werden, 
spricht dafiir, dass diese Systeme physikalisch von hohem Interesse sind. Den Beitrag dieser 
Arbeit dazu will ich hier nocheinmal zusammenfassen: 

• In Kapitel 12 haben wir den Einfluss der Unordung auf die Temperaturabhangigkeit der 
Paarkorrelationsfunktion C{x) = {{ip{x) — (p{Q)Y) im Fall von schwacher {Vipi <C 
hvpcimp) und starker {Vipi » hvpcimp) Unordung untersucht: 

• Im Fall von schwacher Unordnung haben wir mit einem einfachen RG-Argument 
gefunden, dass die Temperaturabhangigkeit dieser Korrelationsfunktion bei kleinen 
Temperaturen konstant ist, d.h. die Unordnung relevant ist, wohingegen bei hohen 
Temperaturen die Korrelationsfunktion linear mit T verlauft, d.h. die Unordnung 
irrelevant wird. Eine quantitativ bessere Erfassung der Temperaturabhangigkeit wurde 
mit Hilfe einer selbstkonsistenten Rechnung im Rahmen einer Burgers-Gleichung erreicht. 
Im Gegensatzzu den Resultaten von Villain und Fernandez [VF84J und Feigel'man [FeiSO] 
haben wir somit ein Verhalten gefunden, welches bei endlichen Temperaturen gilt. 
Zusatzlich wurde das Ergebnis numerisch mit einer Monte-Carlo-Rechnung bestatigt. 

• Im Fall von starker Unordnung konnte die Korrelationsfunktion bei T = 0, unter 
Berucksichtigung der "geordneten" Statistik der Positionen der Verunreinigungen, exakt 
berechnet werden. 

• Der Hauptpunkt der Behandlung des klassischen Modells war die Untersuchung 
der Dynamik. Insbesondere wurde die bisher nicht studierte Kriechdynamik in einer 
Dimension untersucht und folgender Ausdruck fur die Kriechgeschwindigkeit gefunden 



also im Gegensatz zu hoheren Dimensionen ein analytisches Verhalten. Diese Formel 
konnte mit Hilfe von numerischen Simulationen gut bestatigt werden und laBt sich auf 
experimentelle Ergebnisse von S. Zaitsev-Zotov | ,ZZ93j anwenden. 

• Ab Kapitel ^ wurde das klassische Modell zu einem quantenmechanischen mit 
kinetischem Term, welcher in der klassischen Beschreibung nicht berLicksichtigt wurde, 
und einen Term, der Dislokationen bzw. Raum-Zeit-Vortices im Phasenfeld beschreibt, 
verallgemeinert. 

• In den Kapitein |4] und [51 wurden die Phasendiagramme dieses phanomenologischen 
Modells bei tiefen Temperaturen im Fall von eindimensionalen Ladungsdichtewellen 
ausfuhrlich besprochen - mit folgenden Ergebnissen: 

- Bei r = gibt es fiir die Unordnungsfluktuationen einen Lokalisierungs- 
Delokalisierungs-Ubergang bei einem kritischen Wert der Quantenfluktuationen gc^. 
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- AuBerdem findet man fur die Dislokationen zwei Phasen in denen phase-slips 
irrelevant bzw. relevant werden. Der Phasenubergang findet bei einem kritischen 
Wert statt. 1st gc^ldc^ = 24/(pg)^ groBer als 1 hat man drei Phasen: fur kleine 
g eine lokalisierte Phase < ^ci), eine Phase in der sowohl Dislokationen als auch 
Unordnungsfluktuationen irrelevant sind und eine Phase in der phase-slips relevant 
sind {g > gc^). 

— Im Fall der Ladungsdichtewelle ist allerdings gc2 < 5ci- was zur Folge hat, dass 
der Lokalisierungs-Delokalisierungsubergang verschwindet und man direkt von der 
lokalisierten Phase in die phase-slip-bestimmte Phase Obergeht (bei gc2)- 

Die Ergebnisse a us Kapitel |21 wurden teilweise zur Veroffentlichung in Physical Review B 
submittiert. Ein Vordruck dieses Papers ist im Anhang zu finden. 



AbschlieBend will ich kurz darstellen, was ausgehend von den erzielten Resultaten u.a. 
untersucht werden kann: 

• Bisher wurde der Einfluss der phase-slips, wie dies schon in 13.31 angedeutet wurde, 
nicht im Fall eines extern angelegten Feldes untersucht. Wenn phase-slips in dem zu 
betrachtenden System relevant sind, wird dies Konsequenzen fur die Dynamik (z.B. die 
Leitfahigkeit) haben. 

• Man kann zusatzlich noch dissipative EfFekte beriicksichtigen und deren Relevanz fiir das 
System untersuchen. 

• Desweiteren konnen noch Leitfahigkeiten oder Effekte von Tunnelprozessen untersucht 
werden. 
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A Anhang 

A.l Symbole, Einheiten und Notationen 
A. 1.1 verwendete Symbole 



Symbol 


Bedeutung 


Einheit 


Xt = vhf3 


relativistische de Broglie Wellenlange 


L 


Lc 


Fukuya ma-Lee-La nge 


L 


c 


elastische Konstante 




imp 


mittlerer Abstand der Verunreinigungen 


L 


V 


Pinning-Potential-" Dichte" 




X 


Wahrscheinlichkeit zur Ausbildung von 
Dislokationen 


EL-'^ 


kp = Q/2 


Fermi Wellenvektor 


L-i 


A 


(UV-) Cutoff fiir k 


L-i 


{v), VF 


(renormalisierte) Fermi Geschwindigkeit 




P 


Teilchendichte 




K 


Luttinger Parameter 


1 




Phasenfeld und "duales" Feld 


1 


P 


zu (p konjugierter Impulsoperator 




a{x) 


Zufallsphase (gleichwahrscheinlich in [0, 27r[ 
verteilt) 


1 


s = h-^s 


"reduzierte" Wirkung 


1 


n = T-^n 


" reduziertet" Hamiltonian 


1 


i, g, A, V 


dimensionslose Parameter (Bedeutung s. I4.2jl 


1 


Sd = {2TTfKd = 
27r'^/Vr((i/2) 


Oberflache der d-dimensionalen Einheitssphare 






Menge der n x n bzw. n x m Matrizen uber 
einem Korper K 


e(x) 


Stufenfunktion (0 fur x < 0, sonst 1) 



A. 1.2 Notation 

• A: Unordnungsmittelwert 

• (A): thermischer Mittelwert 

• {{A)): Kumulanten-Mittelwert 

• [x]: GauBklammern (Kapitel f2.2.2jl 

• A: Operator / Fouriertransformierte 

• A: Matrix 



A. 2 nutzliche Formein 
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• v;. Vektor {v-w ist das Skalarprodukt oder Matrixmultiplikation) 

• a\b: a modulo b 

• V. imaginare Einheit (um Verwechselungen mit i zu vermeiden) 



Fouriertransformierte und Fourierreihen 

• Die Fouriertransformierte ist in dieser Arbeit wie folgt definiert: 

Sei /(zi) G L^iW^) und r, A; G M^, dann ist die Fouriertransformierte von /(r) gegeben 
durch 

f{k) = [ (frf{r)e-'^-f^ (A.1.1) 

mit der Umkehrtransformation 

f{r) = j^^j^/kf{k)e'^^-^ kurz: jjilDe'-- (A.1.2) 

• Fourierreihen: 

Sei /(r) e L^{L'^) eine periodische /.^-integrierbare Funktion, mit Periodizitatsvolumen 
L'^. Sei k = ^(ni, . . . , n^) mit rij € Z, dann ist die Fourierreihe zu / wie folgt definiert; 

/(i:) = ^"'E^--4 (A-i-3) 

mit J2k = njj=-oo ^""^ "^^"^ Fourierkoeffizienten 

A= / d^rf{r)e'!^-^-. (A.1.4) 

Fiir L ^ oo; 

Afe = 2.L-S (A.1.5) 
d.h. man gelangt wieder zur Fouriertransformierten. 

A. 2 nutzliche Formein 

• Zur Berechnung der n-ten Momente der Abstandsverteilung in Kaoitel l^.2.2l wird 



benutzt. 



• Laurent-Reihe von csch(x) bei x = bis zur 6. Ordnung; 

cschfx = — — ^ = + + Oix') (A.2.2) 

^ ' smh(x) X 6 360 15120 ^ ^ ' ^ ' 
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Zur Berechnung der Korrelationsfunktionen und der RG sind folgende Formein hilfreich: 

OO J 

E ;;2-7;2 = ocotM^I«l) (A-2.3) 



n=— OO ' ' 



/ d^J^ = l_.-\m (A.2.4) 



= . \\ / ^ < X < 27r (A.2.5) 



n=—oo 



V? + a? \a\ sinh(7r|a 
auBerdem fiir die RG des Unordnungsterms 



/•^ 1 

/ dx cosh(ax + 6) = -g [(2 + a^) sinh(a + 6) — 2(a cosh(a + h) + sinh(6))] 

(A.2.6) 



• Poisson'sche Summenformel 



OO OO „ 

E /("^)= E / ^^^/(^)e-^'""" (A.2.7) 

/to 

m=— OO n=— OO 

n-dimensionales GauBintegral: 

/ dVe-5^:^-^-^+^^^ = (27r)'^(detA)-i/2gi2;^.A-i.^;^ ^^2.8) 

wobei A G 7V4„(M) positiv definit und symmetrisch ist und y_ € M". 



Fur ein gauBverteiltes Feld u mit den Mittelwerten [u) = und {u^) ^ folgt mit dem 
Wick-Theorem: 

und daher 

^e^W^ =e(fl^(n))/2_ (A.2.9) 

Insbesondere gilt: 



(2n)! 
2"n! 



2\n 



(sin(p'u)) = 



(A.2.10) 



A. 3 Tilt Symmetrie 
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A. 3 Tilt Symmetrie 



In diesem Anhang wird gezeigt, dass aufgrund der sogenannten statistischen Tilt-Symmetrie 
ip — > ip+ax des Modells ()2.7j) die elastische Konstante hvp nicht renormiert wird. Dazu schaut 
man sich den Replica-Hamiltonian zu ()2.7j) an (siehe dazu Abschnittl4" 



a,/3=l 



1 



d 



dx 



2T 



Irv'o) (1:*=.) (A.3.1) 



y2 ^ 

— cos {ifa - iPf3) } 



In dieser Darstellung sieht man, dass unter o.g. Symmetrietransformation nur ein 
Gradiententerm generiert wird, der durch eine Umdefinition von if wieder "absorbiert" werden 
kann. Nun zeigen wir dies explizit und addieren einen Term —go J^{d(f/dx)dx zu H, wie er 
unter einem Renormierungsgruppenschritt erzeugt werden wLirde, und schreiben den Replica- 
Hamiltonian wie folgt auf: 



dx ( -c 



90-^ COs{ipa - ^p) 



(A.3.2) 



wobei folgende Vektoren und Matrizen verwendet wurden 

A ... i\ 

- 



En 



A = hvpE^ 



V ••• V 

n-dimensionale Einheitsmatrix 



90 = 



{go, . . . ,5o)^ 



Mit Hilfe der Identitaten 

1 , . 1 



a{a + nb) 



[T] , a 7^ —nb und 



2- 
erhalt man 



A. 



Hr. 



--9^LC(n)+ I dx[^-A- 



w 



■ 90^ 



T3^^nd^^2-V2(^ + ^T). 
Die " renormierte" elastische Konstante hvp folgt dann aus der freien Energie: 



wobei C{n) = n/{hvp) + 



1 d^F 



90=0 



(A.3.3) 



(A.3.4) 



66 



Anhang 



Mit {-TZn) = C{n) und F = -Tlim^^o ^{^n - 1) bekommt man schlieBlich 

{hvF)~^ = lim -C(n) = (hvpy^ (A.3.5) 
n^O n 

d.h. die elastische Konstante wird tatsachlich nicht renormiert. 
A. 4 Korrelationsfunktion fiir gauBsches Modell 

In diesem Abschnitt wird die Berechnung der Paarkorrelationsfunktion (((/)(0) — ^(r))^) = 
2 ^^^(0)) — 2 {(f)(0)(p{r)) fiir ein einfaches gauBsches bzw. elastisches Modell 

Ho= f d^rUVrmf (A.4.1) 

kurz beschrieben. Dazu wird zunachst Hq/T in der Zustandssumme um einen linearen 
Quellterm Hq = ^ f d'^rj{r)cl){r) erweitert. Geht man zur Fourierdarstellung (/)(r) = 
L~'^^j^e^--(f)k (entsprechend fiir j) iiber, laBt sich die erweiterte Hamiltonfunktion schreiben 

als {J{k) = ck^) 

H = H,/T + HQ = L~''Y.{^\'t'h?-^hj-^ ■ (A-4.2) 

Mit 4>k = ak+^bk und jk = ak + iPk 'aBt sich die Zustandssumme als Integral iiber die reelen 
Koeffizienten ak und bk darstellen: 

= T\f dak! d6,e-^"1^H+''i)-%"^-''^%) . (A.4.3) 

Nach einer quadratischen Erganzung des Exponeten: a|, — -j^^kf^k = {^o,k ~ 7(^) ~ (^f^) 
(analog fiir 6^) und der Substitution dk = ak — bzw. bk = bk — sind die Integrale 
iiber ak bzw. bk nur noch einfache GauBintegrale, und man erhalt {Zq = Z[j = 0]): 

Z[j] = Zoei^^^J^;^^;^^ . (A.4.4) 

z\j] 

Schreibt man den Exponenten in Z[j] wieder im reaien Raum auf, laBt sich die 
Korrelationsfunktion {(t){r')(f){r)) als Funktionalableitung von Z[j] schreiben: 

(</)(r)0(r')> = j£-Z[j]=TJ-\r-r') 

= TL-<^^e'^(^-^')j-^©. (A.4.5) 
k 

Damit folgt fiir ((^(r') - 0(z:))^) sofort: 



A. 5 Berechnung des Integrals aus 12.2.21 
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Insbesondere fur d = 1 und c = hvp folgt 

((<^(0)-<^(x))2> - 



2r ^ 1 - e'^-^ 
LhvF ^ k'^ 
A 

2T f dk 1 — cos{kx) 



hvp J 27r 

-A 



k^ 



AT 



sm^{kx/2) 



T 



nhvp Jq A^oo hvp 

Beim Ubergang zum Integral wurde ()A.1.5|1 benutzt. 



X . 



(A.4.7) 



A. 5 Berechnung des Integrals aus l2.2.^ 

Berechnung des Integrals: 



Cf= I dxe-'^i— 



2a 



(A.5.1) 



Der quadratische Term in Cfwird zunachst erweitert, was auf die folgenden drei (konvergenten) 
Integra I e fuhrt: 



h = 



j_ r , 2_ 1 



1 

" Jo 



dx e ^x 



2a 



a ^ 



(2fe+l)a 



k=l 



j dx e ^x 


" X ' 




I2a\ 



(A.5.2) 



(2fc-l)o 



(2fe+l)a 



e 



X 



2a 



E 



dx e 



2a 



SO dass C/'= h+h+h- Fur (2fc-l)a < x < (2A;+l)a gilt: [^] A;, d.h. die GauBklammern 
in I2 und konnen durch k ersetzt werden: 

(2k+l)a 



12 = -—'^k J dxe ^x 

13 = ^k'^ I dxe-^ . 



k=l 



(A.5.3) 



(2fc-l)a 



Diese zwei einfachen Integrale liefern 

(2fc+l)a 



y dxe ""x = 2e ^^"((1 + 2A;a) sinh(a) - a cosh(a)) 



(2fc-l)a 
(2fc+l)a 



y dxe"^ = 2e-2^" smh(a) 



(A.5.4) 



(2fc-l)a 
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womit nur noch Summen der Form 

oo oo 

^^6-2^=" und ^fc^e-^^", (A.5.5) 

k=l k=l 

berechnet werden mLissen. Da diese aber nur Ableitungen der geometrischen Reihe sind, sind 
diese auch einfach zu berechnen: 



2ka 



k=l 
oo 



-2a\2 



g_2a _^ g-4a 



k=l 

Man erhalt also letztendlich fur Cf. 



2a \ a sinh(a 

Unter Benutzung von ()A.2.2|) erhalt man folgende Naherungsformel: 

1 7 "^1 
= — - — + ^^a^ + 0(a6) . (A.5.8) 
12 720 30240 ^ ^ ^ ' 

A. 6 Statistik geordneter ZufallsgroBen 

Seien Xi, X2, . . ■ , Xj\[ ZufallsgroBen mit der kumulierten Verteilungsfunktion (cdf) Fi{x) = 
P'^ {Xi < x} wobei im folgenden F{x) = Fi{x) = ... = F/v(x) gelten soil. 

X(2), • • • , ^(AT) seien die geordneten ZufallsgroBen (^(1) < ^(2) < ••• < ^{N)) mit 
der cdf 

n{x) = Pr < ^} = E (I) - n^)f~'' (A.6.1) 

was bedeutet, daB mindestens i der X}^ kleiner oder gleich x sind. In den Spezialfallen i ^ 1,N 
gilt offensichtlich 

p^{x) = Pr{X(;v) < a;} = Pr{alleXi < = 

Pi(x) = Pr{X(i) <x} = l-Pr{alleXi >x} = l-(l-F(x))^ (A.6.2) 

Pi{x) laBt sich auch schreiben als Pi{x) = N — i + l; Pi{x)), wobei I{a, b; x) die normierte 
unvollstandige Beta-Funktion ist, definiert durch 

B{a,b;x) 

I[a,b;x) — 



B{a,b) 

B{a,b-x) = [ f-^l-tf-Ut (A.6.3) 
Jo 

wobei -B(a, b) = B{a, 6; 1) die Beta-Funktion ist. Die Wahrscheinlichkeitsdichte-Funktion (pdf) 
Pi{x) fur kann nun auf zwei verschiedene Arten berechnet werden. Zum einen laBt sich 
die Wahrscheinlichkeit ein X^j) in einem kleinen Intervall [x,x + 6x] zu finden, schreiben als 



A. 6 Statistik geordneter ZufallsgroBen 
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Pi{x)6x = Pi{x + 5x) — Pi{x) + 0{6x'^). fur Sx ^ ist also Pi{x) = Pl{x), womit sofort 
folgt 

^'^(^) = - F{x)f-'f{x) (A.6.4) 

wobei f{x) = F'{x) die pdf der Xj^ ist. Zum anderen kann man sich auch eine mehr 
anschauliche Herleitung fiir pi{x) uberlegen: Wieder soil die Wahrscheinlichkeit fur x < < 
x + 6x bestimmt werden, d.h. die Wahrscheinlichkeit dafur, daB 

• fur i — 1 der Xj. gilt: X^ < x 

• ein Xk € [x,x + 6x] ist und 

• fiir N — i der Xk gilt: X^ > x + Sx 

Fiir diese Aufteilung gibt es offensichtlich (t_i)ijT(jv-i)! =^(jv^i) — B(i N-i+i) Moglichkeiten 
und die Wahrscheinlichkeit fur eine solche Aufteilung ist 

(Pr {Xk < x}y-^ ■Pr{x<Xk<x + 6x}- (Pr {Xk > x + Sx})^'' = 
F^-'^{x) ■ {F{x + 6x) - F{x)) ■ (1 - F{x + 6x))^~' . (A.6.5) 

1 

B{i,N + ' 

F'-^{x) ■ {F{x + 6x) - F{x)) ■ (1 - F{x + fe))^(A.6.6) 
und fiir 6x —>■ folgt wieder der Ausdruck (|A.6.4|) fiir Pi{x). 

Um spater die pdf fiir Differenzen (oder auch anderen Funktionen) von georneten ZufallsgroBen 
zu bestimmen, braucht man noch die gemeinsame pdf fiir zwei GroBen und mit 
X < < X + 6x und y < < y + Sy, wobei 1 < i < j < gelten soil. Diese sei mit 
Pij{x,y) bezeichnet. Geht man analog zu der oben beschriebenen anschaulichen Herleitung 
vor, erhalt man direkt 

P^Jix,v) = _ _ f: _ F'-^-\x)f{x){F{y) - F{x)y-'-'f{y){l - F{y)f-^ 

V ' 

(A.6.7) 

da hier gilt 

• — 1) Xk sind < X 

• ein Xk ist in [x^x + 5x] 

• {j — i — 1) Xk sind in [x + 5x, y] 

• ein Xk ist in [y,y + 5y\ und 

• {N -j) Xk sind > y + 5y 



Ingesamt erhalt man also 

Pr{x < < X + 5x} = 
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d.h. es gibt (j_i)!i!(j_^'i)!i!(jv-j)! Verteilungsmoglichkeiten, jede mit der Wahrscheinlichkeit 
F^-^x){F{x + 5x) - F(x))(F(y) - F{x + Sx)y''-^ {F{y + 5y) - F{y)){l - F{y + 5y)f-^. 
Diese Herleitung kann naturlich auch auf mehr als zwei ZufallsgroBen verallgemeinert werden. 
Jetzt kann die pdf der GroBe = — bestimmt werden. Dazu setzen wir y = Sij+x 
(d.h. Sij ist der Abstand von x und y). Die pdf von ergibt sich sofort aus 

Pi^ij) = / dxpij{x,5ij + x) (A. 6. 8) 

Jut 

Jetzt berechnen wir p{dij) noch in folgendem Spezialfall: f{x) = | '^^^t ' 
d.h. im folgenden seien alle Funktionen nur auf [0,1] von Null verschieden. Fiir die 
Wahrscheinlichkeitsdichte-Funktionen gilt dann 

Pij{x) = Cijx'-\y-xy-'-\l-y)''-^ 



Man erhalt also 



P{Sij) = Cij " dxx'-^6jp-\i-x-6ij)'^-^ (A.6.9) 



Pi^^J) = ^7— n^F^-— r-^C(i - s^Jr-' (A.6.10) 



B{j-i,N-j + i + l) 
und insbesondere im Fall j = i + 1 {6i = 

A. 7 Programm zur Simulation der Dynamik 

Hier ein Listing des Haupteils des Simulationsprogramms: 

//************ program begins . . . ***************** 
//open output-file 

f =f open ( output , " wb " ) ; 

fprintf (f , "N=%d, Saiiiples=°/od, iseed=°/.d\n" ,N,Nsamp, iseeda) ; 
fprintf (f , "time=y.le, dt=y.le, T=y.le\n" ,time,dt,T) ; 

f or(k=0;k<nE;k++) // loop over driving force 
{if(nE==l) E=(bE+vE)/2; 
else E=vE+k*(bE-vE)/(nE-l) ; // electric field or driving force 



for(j=0; j<Nsamp; j++) 
{iseeds=iseeda+j *108; 
x=ran2(iseeds) ; //only for initialzation 
//Initialization 
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for(n=0;n<N;n++) 
{phi [n] =0 . ; //initially surface is flat 

// random phase: a[n] 
a[n] =M_2PI*ran2() ; 

} 

Jsav=0.0; // to calculate average J_cdw 
avc=0; //counter for " 
//TIME EVOLUTION 
f or (i=0; Knstep; i++) 
{f or(n=0;n<N;n++) 

{ np=n+l;if (np==N) np=0; 
nm=n- 1 ; if (nm<0 ) nm=N- 1 ; 

noise [n] =Tsdev*sqrt (-2 . 0*log(ran2 () ) ) *cos (M_2PI*ran2 () ) ; 

//calc d(phi)/dt 

deril [n] =(phi [np] +phi [nm] -2*phi [n] )+sin(phi [n] -a[n] ) 
//discrete laplacian pinning 

-E +noise [n] ; 

//driving force temperature 

} //end n 

f or(n=0;n<N;n++) phit [n] =phi [n] +dt*deril [n] ; 

//modified RK: 
for(n=0;n<N;n++) 
{ np=n+l ; if (np==N) np=0; 
nm=n-l ; if (nm<0) nm=N-l; 

deri2 [n] =(phit [np] +phit [nm] -2*phit [n] )+sin(phit [n] -a[n] )-E+noise [n] ; 
} //end n 
f or(n=0;n<N;n++) 
{phi [n] = (phi [n] +phit [n] ) /2 ; 
deri [n] = (deri 1 [n] +deri2 [n] ) /2 ; 

} 

for (n=0 ; n<N ; n++) phi [n] =phi [n] +dt*deri [n] ; 

if ( (i>=nskip) && (i°/oical==0) ) // calculate the average J_cdw 
{Jav=0.0; 
f or(n=0;n<N;n++) Jav+=deri [n] ; 
Jav=Jav/N; 
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Jsav+=Jav; 
avc++; 

> 

}//end time 
Jsav= Jsav/ avc ; 
savj [j] =Jsav; 
}//eiid sample-loop 
// sample average for J_cdw 
J 1=0.0; // final average 
J2=0.0; 

f or(j=0; j<Nsamp; j++) 

{Jl+=savj [j] ; J2+=savj [j] *savj [j] ; } 
// calculate error bars due to sample average 
if (Nsamp>l) err=sqrt ( ( J2- Jl* Jl/Nsamp) / (Nsamp* (Nsamp-1) ) ) ; 
else err=0.0; 

fprintf (f , "7,le "/.le y,le\n" ,E, Jl/Nsamp, err) ; 
} // end of loop over driving force 
f close (f ) ; 
//finish 

Bemerkungen: 

• Die Funktion ran2() liefert eine Zufallszahl im Intervall [0,1] (entnommen 
aus (EIVE92J) 

• Tsdev=sqrt (2 . 0*T) ; ist die Standardabweichung fur die gauBverteilten Zufallszahlen 

viT) 
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A. 8 Teilpublikation 

Andreas Glatz and Mai Suan Li, Collective dynamics of one-dimensional charge density waves, 
Phys. Rev. B 64, 184301 (2001). (Teile von Kapitel 2) 

Andreas Glatz and Thomas Nattermann, One-Dimensional Disordered Density Waves and 
Superfluids: The Role of Quantum Phase Slips and Thermal Fluctuations, Phys. Rev. Lett. 88, 
256401 (2002). (basierend auf Kapitel 3 und 4) 
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